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ABSTRACT 


The  behavior  of  steady  flow  in  two  successive  90°  pipe 
elbows  of  large  dimensionless  curvature  (  a/R  =  0.357  )  has 
been  investigated  experimentally.  The  flow  has  been  establi¬ 
shed  in  terms  of  the  dimensionless  parameters;  the  Reynolds 
(or  Dean.  )  number,  Re  or  D,  the  non-dimensional  length  of  the 
spacer  between  the  elbows,  1',  and  the  angle  between  the  el¬ 
bows'  planes  of  curvature,  V. 

The  influence  of  the  dimensionless  variables  on  the  la¬ 
minar  secondary  flow  has  been  shown  for  the  case  of  uniform 
velocity  profile  at  the  entrance  and  the  types  of  the  secon¬ 
dary  flow  patterns  have  been  classified.  Also  the  association 
of  the  flow  types  to  the  two  basic  flow  structures,  observed 
when  the  elbows  have  formed  the  U  shape  and  the  S  shape,  has 
been  given. 

The  positions  of  the  separation  lines  and  the  reattach¬ 
ment  lines  have  been  measured  and  presented  as  functions  of 
the  dimensionless  variables,  both  for  the  uniform  velocity 
profile  and  the  fully  developed  flow  at  the  entrance. 

The  flow  inside  the  separation  zones  (backflow),  the 
influence  of  the  second  elbow  on  the  flow  upstream  and  the 
transition  to  turbulence  also  have  been  observed. 
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CHAPTER  1 


INTRODUCTION 

1.1  Introductory  Remarks 

Problems  associated  with  fluid  motion  through  two  or 
more  elbows,  have  attracted  a  good  deal  of  attention  re¬ 
cently,  mainly  because  of  their  possible  relevance  to 
many  engineering  or  physiological  situations.  For  example, 
a  pipe  line  in  practical  use  often  has  many  bent  portions 
located  close  to  one  another  and  the  flows  through  these 
elbows  are  affected  by  their  relative  positions.  In  this 
case,  the  total  bend  loss  is  not  equal  to  the  sum  of  the 
individual  elbow  losses  as  would  be  the  case  if  the  elbows 
were  located  separately  sufficiently  far  apart  in  a  pipe 
line.  This  is  due  to  the  mutual  interference  of  the 
elbows  on  the  flow. 

Another  possible  application  of  this  work  is  the 
blood  flow  in  the  cardiovascular  system:  two  90° elbows 
with  no  spacer  and  angle  ¥  =  15° can  simulate  the  human 
aortic  arch.  In  contrast  to  the  interests  of  engineering, 
where  the  increase  in  resistance  due  to  curvature  is  of 
principal  interest,  a  knowledge  of  the  velocity  distribution 
is  required  for  the  study  of  the  cardiovascular  system. 
Through  a  knowledge  of  the  velocity  distribution,  the 
distribution  of  injected  substances  and  the  effects  of  the 
distribution  of  wall  shearing  stress  on  the  formation  of 
arterial  lesions  may  be  better  understood. 
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1.2  Literature  Survey 


Only  a  few  studies  on  the  subject  of  fluid  flow  in 
pipes  with  two  or  more  bends  combined  have  been  carried  out, 
but  much  attention  has  been  directed  to  the  more  general, 
and  theoretically  simpler  problem  associated  with  the  fluid 
motion  through  a  curved  tube. 

The  first  investigation  of  loss  of  head  in  a  coiled  tube 
was  done  by  Eustice  ( 1 9 1 1 ) ,  who  carried  out  color  band  ex¬ 
periments  and  classified  the  flows  as  laminar,  double  hel¬ 
ical  and  turbulent.  Eecause  of  the  elliptical  cross-section 
of  the  tube  used  (due  to  the  bending  process  employed)  the 
results  obtained  for  head  loss  were  of  little  practical  sig¬ 
nificance  . 

The  theoretical  behavior  of  the  steady  motion  of  a  vis¬ 
cous  fluid  in  a  curved  tube  having  a  circular  cross-section, 
was  first  considered  by  Dean  (1927 t  1928).  In  this  analysis 
it  was  found  necessary  to  assume  that  the  dimensionless  cur¬ 
vature  of  the  pipe  (ratio  a/R)  is  small.  In  the  first  ap¬ 
proximation  (by  Fourier  -  series  analysis  method  )  to  the 
Navier-Stokes  equations  (Dean,  1927)  the  relation  between 
the  rate  of  flow  and  the  pressure  gradient  was  found  to  be 
independent  of  the  curvature  of  the  tube,  but  in  a  higher 
order  analysis,  (Dean  1928,  fourth  approximation)  the  re¬ 
duced  flow  due  to  curvature  was  determined  to  be  a  function 
of  a  single  constant  K:=|-Re2a/R.  The  problem  of  determining 
the  relation  between  flow  rate  and  pressure  gradient  in 
any  curved  pipe  of  small  dimensionless  curvature  therefore 


. 


reduced  to  that  of  finding  a  single  function  f (K) .  Dean 
showed  that  his  analysis  was  in  good  qualitative  agreement 
with  Eustice's  (1911)  experiment,  reliable  for  values  of  K 
less  than  576, and  gave  an  equatiori  for  the  axial  pressure 
drop.  This  analysis  cannot  be  used  to  determine  the  vel¬ 
ocity  distribution  which  is  important  in  many  engineering 
problems . 

White  (1929)  continued  the  experimental  investigation 
of  laminar  motion  in  a  curved  pipe  for  fully  developed  en¬ 
trance  profile  and  small  dimensionless  curvature.  In 
his  paper  an  empirical  equation  for  the  axial  pressure  drop 
in  a  coiled  pipe  was  suggested.  Also,  limitations  on  the 
application  of  Dean's  theory  and  a  comparison  with  Eustice' 
experiment  were  given.  White  concluded,  that  even  with  the 
small  dimensionless  curvature  of  a  pipe  (a/R=l/50)  flow  was 
fully  laminar  up  to  Re=6000.  This  result  was  experiment¬ 
ally  verified  by  Taylor  (1929)  by  repetition  of  the  Eustice 
experiment,  but  with  focus  on  turbulence  (three  pipes  of 
small  dimensionless  curvature  were  used) . 

Adler  (1934)  showed  that  the  ratio  of  the  resistance 
coefficients  for  the  curved  and  a  straight  pipe  of  the  same 
radius,  depends  only  on  K  (the  constant  introduced  by  Dean) 
as  long  as  the  motion  is  laminar. 

Based  on  the  White,  Taylor,  and  Adler  experiments,  se¬ 
veral  attempts  have  been  made  to  extend  the  solution  over 
the  entire  range  of  laminar  flow. 

Barua  (1963)  deduced  from  experimental  results  that 
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for  large  K  , motion  outside  the  boundary  layer  is  mostly 
confined  to  planes  parallel  to  the  plane  of  symmetry  of  the 
pipe.  This  concept  led  to  the  prediction  that  the  effect 
of  curvature  is  to  increase  the  resistance  coefficient  of 
the  curved  pipe  relative  to  that  of  the  same  pipe  if  it  was 
straight,  according  to  the  empirical  formula 

res.  coeff_. _ curved  pipe _ =  0>0918  Ki  +  0>509 

res.  coeff.  straight  pipe 

Barua's  solution  was  obtained  using  the  momentum  integral 
method . 

Topakoglu  (1967)  gave  a  systematic  method  for  determin¬ 
ing  an  approximate  solution  (with  required  accuracy)  for 
steady  laminar  flow  in  a  curved  pipe.  Using  the  assumption 
that  the  stream-function  and  the  normal  component  of  the  vel¬ 
ocity  vector  may  be  expanded  in  a  power  series  in  terms  of 
the  curvature  of  the  pipe  he  obtained  an  expression  for  the 
rate  of  flow  in  a  curved  pipe. 

McConalgue  &  Srivastava  (1968)  extended  Dean's  work 

JL 

and  adopted  a  parameter  D=2Re ( 2a/R) 2 .  Solutions  obtained 
numerically  by  Fourier  analysis,  were  given  over  the  range 
D=96  to  605.72,  the  value  D=96  corresponding  to  the  upper 
limit  K=576  of  Dean's  work.  Truesdell  and  Adler  (1970)  ob¬ 
tained  results  up  to  D=35?8  by  a  method  based  on  finite-dif¬ 
ference  approximations  to  the  governing  partial  differential 
equations.  By  a  similar  method,  Greenspan  (1973)  obtained  a 
numerical  solution  covering  the  whole  range  of  D  for  laminar 
flow,  that  is,  for  0<D^  5000.  In  his  paper,  calculations 
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applied  to  the  particular  model  studied  by  Dean  and 
McConalgue  &  Srivastava  and  velocity  profiles  for  different 
D  were  given.  However,  as  was  pointed  out  by  Yao  &  Berger 
(1975)  the  numerical  viscosity  inherent  in  the  finite-differ¬ 
ence  scheme  may  be  different  than  the  actual  viscosity  and 
therefore  the  results  may  not  be  valid  for  large  D. 

Earlier,  Ito  (1969),  based  on  assumptions  similar  to 
those  used  by  Barua  (1963)*  investigated  theoretically  the 
laminar  flow  in  a  curved  pipe  and  derived  a  formula  for  the 
friction  factor.  Also  a  comparison  of  his  calculated  results 
with  the  experimental  results  of  Adler  (193^)  and  Ito  (19 59 ) 
was  given. 

By  the  method  of  successive  approximations,  Sankaraiah  & 
Rao  (1973)  solved  the  equation  of  flow  obtained  in  terms  of 
a  secondary  flow  stream  function  and  an  axial  velocity  com¬ 
ponent  (as  suggested  by  Dean).  The  first  five  approximations 
were  computed  and  the  solution  obtained  was  used  to  deter¬ 
mine  the  axial  velocity  distribution,  secondary  flow  patterns, 
axial  pressure  drop  and  pressure  distribution  along  the  pipe 
wall.  A  semiemprical  equation  was  obtained  for  axial  pres¬ 
sure  drop  and  the  theoretical  results  were  compared  with  the 
available  experimental  data. 

Collins  &  Dennis  (1975)  attempted  to  establish  an  asymp¬ 
totic  structure  for  the  behavior  of  the  numerical  solutions 
as  d  -roo.  The  solutions  were  obtained  by  the  method  of  for¬ 
ward  and  backward  differences  and  corrected.  The  complete 
range  covered  by  McConalgue  &  Srivastava  (1968)  and 
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Greenspan  (1973)  was  considered  and  the  results  were  found 
to  be  generally  consistent  with  the  former  but  quantitative¬ 
ly  different  from  the  latter,  although  many  of  the  general 
trends  were  qualitatively  similar*  The  results  obtained 
were  in  good  agreement  with  their  experiments. 

Smith  (1976)  investigated  the  influence  of  curvature  on- 
fluid  flow  for  a  pipe  that  starts  bending  uniformly  after  an 
initial  straight  section.  The  Reynolds  number  and  dimen¬ 
sionless  curvature  were  assumed  to  be  large  and  small,  re¬ 
spectively.  Smith's  theory  was  also  extended  to  larger 
length  scale,  larger  D  and  to  various  velocity  profiles. 

Patankar  et  al  (197^)  applied  a  calculation  proce¬ 
dure  (finite-differences  marching  technique)  to  three-dim¬ 
ensional  parabolic  flows  to  predict  the  velocity  and  temp¬ 
erature  fields  in  helically  coiled  pipes,  in  the  developing 
and  fully  developed  regions.  Results  were  compared  with  ex¬ 
perimental  data  and  showed  good  agreement.  Also,  the  ef¬ 
fects  of  Dean's  number  on  the  friction  factor  were  present¬ 
ed.  Only  the  small  dimensionless  curvatures  (  a/R  <.0.07  ) 
were  considered. 

Murata  et  al  (1976)  investigated  theoretically  the 
laminar  flow  through  pipes  of  circular  cross-section,  but 
with  locally  varying  curvature  of  the  centre  line.  The 
analysis  is  applicable  to  any  two-dimensional  curved  pipe, 
when  the  centre-line  dimensionless  curvature  is  small. 

Time  dependent  fully  developed  flow  in  curved  pipes  was 
analysed  by  Lyne  (1970),  Zalosh  &  Nelson  (1973) »  Smith  (1975) 
and  Singh  et  al  (1978). 
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All  of  the  theoretical  investigations  mentioned  above 
were  carried  out  with  the  assumption  of  small  dimensionless 
curvature  of  a  pipe  (0<  a/R<l/50).  This  assumption  allow- 

?  O 

ed  the  investigators  to  neglect  terms  of  order  a  /R  and 
significantly  simplified  the  calculations,  but  for  practical 
applications,  even  for  pipes  of  small  dimensionless  curva¬ 
ture,  results  obtained  in  this  way  were  sometimes  unsatis¬ 
factory  (see  McConalgue  &  Srivastava  1968,  Greenspan  1973 
and  Collins  &  Dennis  1975) •  Although  any  solution  obtained 
with  the  above  assumption  would  not  be  valid  for  bends  of 
large  dimensionless  curvature,  some  investigators  have  ap¬ 
plied  it  to  an  entry  region  of  curved  pipe. 

By  the  method  of  matched  asymptotic  expansion,  Singh 
(197*0  obtained  the  solution  for  flow  development  in  a  curved 
tube  near  the  inlet  for  conditions  of  constant  dynamic 
pressure  at  the  entrance  (application-blood  flow  in  an  aorta) 
and  uniform  velocity  profile.  This  solution,  however, 

i 

breaks  down  at  a  distance  of  the  order  (aR)2,  corresponding 
physically  to  the  point  beyond  which  the  effect  of  the  cen¬ 
trifugal  force,  initially  small,  becomes  as  important  as 
inertia  and  viscous  forces  within  the  boundary  layer. 

Yao  &  Berger  (1975)  investigated  flow  with  large  D  as 
it  developed  from  a  uniformly  distributed  velocity  at  the 
entrance  to  a  fully  developed  profile.  Barua's  results  for 
the  fully  developed  flow  were  adopted  as  down-stream  condi¬ 
tions  . 

These  theoretical  considerations  were  inspected  experi- 
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mentally  by  Agrawal  et  al  (1978)  applying  laser  anemometry 
for  the  measurements  of  the  axial  velocity  and  the  component 
of  the  secondary  velocity  (parallel  to  the  plane  of  curva¬ 
ture)  .  They  dealt  with  the  laminar  flow  of  a  Newtonian 
fluid  in  the  entry  region  of  a  curved  pipe  of  curvatures 
1/20  to  l/7,  covering  Dean's  numbers  from  138  to  679.  The 
experimental  axial  velocity  profiles  were  compared  with 
those  constructed  from  the  theoretical  analysis  of  Singh 
(197*0  and  Yao  &  Eerger  (1975)*  The  quantitative  agreement 
between  theory  and  experiment  was  found  to  be  poor.  However, 
some  of  the  features  observed  in  the  experiment  were  in  qua¬ 
litative  agreement  with  the  theoretical  solution  of  Yao  & 
Berger . 

Singh  (1978),  continuing  his  study  of  the  problem  with 
particular  reference  to  the  ascending  aorta,  investigated 
the  oscillatory  flow  of  viscous  fluid  in  the  entrance  to  the 
curved  pipe.  He  neglected  arterial  wall  elasticity,  the 
taper  of  the  aorta  and  the  influence  of  the  branches. 

A  different  approach,  more  important  from  the  point  of 
view  of  this  work  because  of  possible  application  to  the 
large  curvature  elbows,  was  presented  by  another  group  of 
investigators.  They  carried  out  a  theoretical  analysis  of  an 
inviscid  flow  without  body  forces,  which  significantly  re¬ 
duced  the  mathematical  difficulties.  First,  Squire  &  Winter 
(1949,  1951)  showed  theoretically  that  secondary  flows  could 
occur  in  a  bend  through  which  a  perfect  (inviscid)  fluid  is 
flowing  as  a  result  of  a  non-uniform  distribution  of  velocity 


at  the  entrance  to  the  bend. 

Hawthorne  (1950)  applied  Squire  &  Winter's  conclusion  to 
a  more  general  theoretical  investigations  of  the  rotational 
flow.  It  was  found  that  the  secondary  flow  was  not  spiral  hut 
oscillatory, and  that  the  direction  of  the  circulation  chang¬ 
ed  periodically.  Hawthorne  showed  also  that  the  secondary 
circulation  remains  unchanged  along  the  streamline  if  the 
direction  of  acceleration  (or  pressure  gradient)  lies  in  the 
plane  containing  the  velocity  vector  and  the  normal  to  the 
Bernoulli  surface. 

Later,  Horlock  (1953)  described  some  experiments  on  the 
secondary  flows  in  curved  pipes,  dealing  in  particular  with 
the  effects  of  reversal  of  curvature  on  the  pattern  of  the 
streamlines . 

Recently,  Murakami  &  Shimizu  (1978)  extended  Horlock’ s 
analysis  for  three  large  dimensionless  curvature  elbows  com¬ 
bined  in  series  and  compared  the  results  with  experiments. 

The  investigators  predicted  three  types  of  swirling  flows  to 
occur:  oscillatory,  spiral  and  suspended.  Experiments  were 

carried  out  for  fully  developed  turbulent  profile  (Re=l0^)f 
different  H*  angles  and  different  lenths  of  the  spacers  be¬ 
tween  the  elbows. 

Due  to  the  difficulty  of  obtaining  even  approximate  theo¬ 
retical  solutions,  some  investigators  have  focused  their  atten¬ 
tion,  for  practical  purposes,  on  purely  empirical  work.  Those 
included  are  Beij  (1938) »  Crane  Ltd.  (19^2),  Pigott  (1950) » 

Ito  ( 1959 »  I960)  and  others. 
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The  hydraulic  losses  and  flow  patterns  in  pipes  with  two 
bends  combined  for  fully  developed  turbulent  profile  were 
investigated  by  Murakami  &  Shimizu  (1977). 

1.3.  Objectives  of  the  Thesis 

The  aim  of  this  thesis  is  to  provide  a  better  under¬ 
standing  of  the  flow  of  fluid  through  the  90°  elbows  locat¬ 
ed  close  to  one  another.  The  variable  parameters  will  be: 

-  the  Dean  number 

-  the  length  of  the  spacer  between  the  elbows 

-  the  angle  between  the  planes  of  curvature  of  the  elbows 

Qualitative  and  quantitative  results  of  the  separation 
and  reattachment  positions  will  be  produced  for  two  cases: 
a  constant  pressure  tank  or  a  fully  developed  flow  at  the 
inlet.  Also  qualitative  results  of  the  observations  of  the 
turbulence  and  secondary  flow  patterns  will  be  given  and  il¬ 
lustrated  by  photographs. 
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CHAPTER  II 


DIMENSIONLESS  PARAMETERS 
2.1  Introductory  Remarks 

As  has  been  mentioned,  a  theoretical  approach  presents 
very  serious  computational  difficulties.  Moreover,  there 
are  no  effective  simplifying  assumptions  useful  in  the  case 
being  considered.  Neglecting  the  fluid  viscosity  (e.g. 
Horlock  195d,  Murakami  &  Shimizu  1978),  is  acceptable  for 
high  Reynolds  number  flows  and/or  non-uniform  velocity  pro¬ 
file  flow,  but  for  the  case  being  considered,  this  simpli¬ 
fication  leads  to  a  direct  contradiction  with  the  existence 
of  secondary  flow.  Similarly,  the  assumption  that  viscous 
forces  exist  in  the  inlet  pipe  only  (for  a  non-uniform  vel¬ 
ocity  distribution  at  the  inlet  to  the  first  elbow)  is  also 
unsatisfactory  because  the  inlet  pipe  is  relatively  short 
and  therefore  the  flow  development  is  not  completed  before 
the  inlet  to  the  first  elbow,  even  for  the  lowest  Reynolds 
numbers .  But  the  rapid  increase  of  the  boundary  layer  thick¬ 
ness  in  the  inlet  pipe  also  indicates  the  continuation  of 
this  process  inside  the  elbow  and  the  importance  of  the  vis¬ 
cosity  forces. 

Additionally,  the  large  dimensionless  curvature  of  the 

p 

elbows  means  that  terms  of  the  orders  (a/R)  and  (a/R)  are 
equally  important  and  the  methods  used  by  Dean  (1927,  1928), 
McConal.gue  &  Srivastava  (I968),  Greenspan  (1968)  and  others 
may  not  be  adopted  for  this  study. 

For  the  reasons  described  above  this  thesis  will  focus 
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on  the  expirimental  study  of  the  water  flow  in  two  com¬ 
posite  pipe  elbows  with  special  application  to  the  blood 
flow  in  the  human  aortic  arch,  assuming  that  the  blood  in 
large  arteries  behaves  like  a  Newtonian  fluid. 

2.2  Dimensionless  parameters 

The  purpose  of  this  section  is  to  establish  a  set  of 
dimensionless  parameters  fully  describing  the  system  and  con¬ 
cerned  with  the  conditions  under  which  flows  of  different 
fluids  (e.g.  water  and  blood)  in  two  geometrically  similar 
elbow  systems  are  dynamically  similar. 

The  physical  quantities  which  determine  the  geometry 
and  the  flow  in  the  two  elbows,  are: 

i.  geometry:  radius  of  curvature  R  and  internal  diameter  of 
the  elbows  d;  mutual  positions  of  the  elbows  (length  of 
the  spacer  L  and  angle  between  the  curvature  planes  of  the 
elbows  V  ) 

ii.  flow:  mean  velocity  U  ;  properties  of  the  fluid  (density 
3>  and  viscosity^)  ;  the  representative  linear  dimension 

of  the  system  (diameter  of  pipe  d) 

It  can  be  proved  by  the  condition  of  dynamic  similarity 
or  by  the  method  of  indices  (Schlichting  1968)  that  there 
exists  a  unique  dimensionless  combination  of  the  four  quan¬ 
tities:  ,  <?  and  d  -  the  Reynolds  number  Re,  which  pro¬ 

vides  the  conditions  under  which  flows  of  different  fluids 
about  two  geometrically  similar  bodies  (in  this  case-  in  two 
geometrically  similar  elbow  systems) ,  and  with  identical  in- 
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itial  flow  direction  display  geometrically  similar  stream¬ 
lines  (are  dynamically  similar),  when  only  frictional  and  in¬ 
ertia  forces  are  important. 

For  geometrical  similarity,  the  linear  parameters 
must  also  he  expressed  in  a  dimensionless  form.  This  can  he 
done  hy  referring  them  to  the  representative  linear  dimen¬ 
sion  of  the  system,  the  diameter  d.  Introducing  also 
L'  =  L/d  and  taking  into  account  that  the  angle  ^  is  dimen¬ 
sionless  hy  definition,  the  following  set  of  the  dimension¬ 
less  parameters,  fully  describing  the  system,  have  been  ob¬ 
tained: 

R/d, L' ,  V,Re 

As  has  been  mentioned  in  Section  2.1,  Dean  (1928) 
showed  that  the  dynamic  similarity  of  the  fully  developed 
flow  in  a  single  bend  of  small  dimensionless  curvature  de¬ 
pends  on  a  non-dimensional  parameter,  called  Dean's  num¬ 
ber  by  McConalque  &  Srivastava  ,  which  is  a  function  of  the 
Reynolds  number  and  the  geometrical  parameters  of  the  bend. 
Although  in  this  case  the  flow  is  not  fully  developed  and 
the  dimensionless  curvature  of  the  elbows  is  large,  both 
Reynolds  and  Dean  numbers  will  be  used  in  present  work,  for 
convenience  of  the  reader.  It  should  be  noted,  that  because 
of  the  fixed  geometry  used  in  this  investigation,  the  Dean 
number  is  the  Reynolds  number  multiplied  by  a  constant. 

It  must  be  mentioned  also  that  there  exist  at  least 
three  definitions  of  the  Dean  number.  They  are  as  follows: 
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McConalgue  &  Srivastava  (1968),  Greenspan  (1973) 

D  =  2  Re  (2a/R)^  (2.1) 

Singh  (1974),  Yao  &  Berger  (1975) 

D  =  Re  (a/R)^  .  (2.2) 

Schlichting  (1968) 

D  =  \  Re  ( a/R ) *  (2.3) 

The  definition  (2.1)  which  is  the  most  widely  used,  will  be 
employed  in  this  thesis. 

Finally,  the  set  of  the  dimensionless  variables 
includes : 

-  Reynolds  (Dean)  number, Re  (D) 

-  Dimensionless  spacer  length,  L' 

-  Angle  between  the  curvature  planes  of  the  elbows,  V 


CHAPTER  III 


APPARATUS 

3.1  Apparatus  Description 

The  purpose  of  the  experimental  model  was  to  enable  the 
study  of  the  flow  characteristics  in  two  successive  pipe  el¬ 
bows  and  relate  it  to  blood  flow  in  the  aortic  arch.  For 
this  reason,  the  dimensions  of  the  elbows  and  the  lengths  of 
the  inlet  and  outlet  pipes,  were  chosen  to  be  close  to  the 
average  values  obtained  from  measurements  made  on  humans  of 
the  ascending  aorta,  the  aortic  arch,  and  the  section  of  the 
aorta  downstream  from  this  arch  (McDonald  i960,  Singh  et  al 
1978).  The  aorta  is  an  elastic  tube  that  stems  from  the 
left  ventricle  of  the  heart,  curving  in  a  complicated  three- 
dimensional  way  with  branchings  to  the  head  and  upper  limbs. 
In  this  investigation,  the  aorta  has  been  modelled  using  two 
elbows  with  a  straight  pipe  at  the  inlet  and  a  straight  pipe 
at  the  outlet.  As  a  first  approach  to  the  problem  the  three 
dimensional  model  employed  ignores  the  influence  of  branches, 
the  taper  of  the  aorta,  the  arterial  wall  elasticity  and  the 
pulsating  character  of  the  blood  flow.  The  various  length 
spacers  have  been  introduced  between  the  elbows  to  make  the 
model  more  universal  and  the  results  obtained  more  general 
and  useful. 

As  can  be  seen  in  Figure  1 ,  some  additional  devices 
have  been  necessary  to  provide  the  working  conditions  for 
the  model.  To  obtain  the  uniform  velocity  profile  at  the  in¬ 
let  to  the  first  elbow,  city  water  was  supplied  to  the  clo- 
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sed,  watertight  inlet  tank.  For  the  uniform  velocity  pro¬ 
file,  the  long  copper  tube  was  used,  but  due  to  supply  dif¬ 
ficulties  the  internal  diameter  of  this  tube  d^  differed 
slightly  from  the  diameter  of  the  elbows  (see  Table  1)  and 
an  adapter  had  to  be  used.  Therefore  the  velocity  profile 
produced  during  these  tests  cannot  be  considered  to  be  fully 
parabolic,  but  as  being  very  close  to  it.  The  required  length 
of  this  pipe  was  calculated  on  the  basis  of  the  formula, 
(Daugherty  &  Ingersoll  195*0 

L  p  =  0.058  d  p  Re  (3-D 

and  provided  fully  developed  laminar  flow  at  the  outlet  of 
the  pipe  up  to  Re  =  1918. 

To  avoid  thermal  effects  the  water  temperature  was 
maintained  at  the  same  level  as  the  room  temperature  by  mix¬ 
ing  hot  and  cold  water.  The  flow  rate  was  controlled  and 
the  safe,  Tow  pressure  in  the  supplying  pipe  was  maintained 
using  the  flow  control  valves.  The  flow  rate  was  coarsely 
measured  by  the  flow  meter  and  precisely  using  the 
scaled  flow  rate  measurement  tank  and  a  stopwatch. 

The  outlet  unit,  containing  the  outlet  tank,  the  flow 
rate  measurement  tank  and  the  bottom  tank,  were  mounted  on 
casters,  to  enable  easy  changes  of  the  angle  ^ .  This  angle 
was  changed  by  loosening  the  vertical  screws  pressing  the 
elbows  together,  located  in  the  steel  housing  of  the  elbows, 
and  rotating  the  unit  about  the  vertical  axis  of  the  spacer. 
Additionally,  the  outlet  tank  located  at  the  unit  could  be 
moved  vertically  to  keep  the  outlet  pipe  in  the  horizontal 


position,  when  the  spacer  length  was  changed. 

Apparatus  data  are  provided  in  Table  1.  The  general 
views  and  the  details  of  the  apparatus  are  shown  in  Figures 
2  through  5  inclusive. 

3.2  Flow  Visualisation 

The  flow  was  made  visible  by  injecting  water  colored 
with  a  dark-blue  dye  using  dye  pipes  located  at  various 
points  from  the  inlet  tank  to  the  outlet  pipe.  For  the 
exact  dye  pipe  locations  see  Table  2.  The  first  dye  pipe, 
marked  A,  was  mounted  (see  Figures  1,  3,  &  4)  such  that  it 
could  be  located  at  any  point  before  the  entrance  to  the  in¬ 
let  pipe  in  the  inlet  tank.  That  is,  it  was  possible  to  in¬ 
troduce  the  dye  mixture  at  any  point  of  the  entrance.  All 
other  dye  pipes  were  fixed  perpendicularly  to  the  main  di¬ 
rection  of  flow  in  the  walls  of  the  elbows,  spacer  and  the 
outlet  pipe  (see  Figures  2  &  4)  .  Due  to  the  fixed  position 
of  the  dye  pipes  located  in  the  elbows  walls,  two,  three,  or 
four  dye  pipes  were  located  at  the  same  cross-section  and 
were  marked  left,  right,  outer  and  inner. 

Dye  pipes  being  in  use  during  a  certain  experiment  were 
connected  to  the  container  supplying  the  dye  mixture  and  all 
other  dye  pipes  were  closed. 

The  dye  motion  in  the  elbows  was  observed  and  photo¬ 
graphed,  and  to  obtain  the  true  three-dimensional  pattern 
of  the  flow,  two  small  variable  position  mirrors  were  used 
(see  Figure  6).  The  vertical  position  of  the  mirrors  was 
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in  accordance  with  the  needs  of  the  experiment  being  con¬ 
ducted.  Since  the  elbows  were  made  by  matching  the  two 
halves  then  joining  them  together  with  the  connecting  plane 
coinciding  with  the  plane  of  curvature  of  the  elbows,  care 
was  taken  to  position  the  mirrors  to  eliminate  possible  op¬ 
tical  effects  created  by  this  joint,  see  Figure  6. 

3.3  Apparatus  Limitations 

The  actual  dimensions  of  the  elbows  chosen  here  intro¬ 
duced  certain  limitations  to  the  experimental  procedure  which 
will  be  described  in  Chapter  IV.  This  is  related  to  the 
relatively  small  diameter,  from  the  experimental  point  of 
view,  that  has  prevented  the  locating  inside  the  elbows  or 
pipes  of  any  pitot  tubes  or  other  velocity  measuring  devices 
of  the  diameter  larger  than  approximately  1  mm.  This  means 
that  for  low  velocity  flows  the  dynamic  force  driving  a  pres¬ 
sure  sensitive  device  connected  to  a  pitot  tube  will  be  very 
low.  For  example,  for  the  more  advantageous  case  of  poten¬ 
tial  water  flow  with  a  velocity  of  4.5  cm/s,  the  average  vel¬ 
ocity  for  Re  =  500  (D  =  2535)  in  this  experiment,  the  dyna¬ 
mic  force  on  the  1  mm  diameter  pitot  tube  and  the  dynamic 

pressure  head,  using  the  Bernoulli  equation,  would  be  ap- 

-6) 

proximately  0.8  x  10  N  and  0.1  mm  of  water,  respectively. 
These  numbers  are  lower  for  viscous  flow  and  decrease  rapid¬ 
ly  as  the  flow  velocity  decreases  (0.005  mm  H^O  for  1  cm/s). 
To  investigate  secondary  flow  velocities,  which  are  at  least 
one  order  of  magnitude  smaller  than  the  main  flow  velocities, 


the  necessary  sensitivity  of  the  device  measuring  the 
differential  pressure  must  he  better  than  0.001  mm 
H^O.  No  device  has  been  found  to  satisfy  this  condition  - 
the  very  sensitive  condenser-type  eensors,  designed  to  work 
with  gases,  require  a  special  interfacing  device  for  use  with 
water,  which  reduces  the  sensitivity  to  the  point  where  it 
is  of  no  use . 

For  the  reasons  noted  above,  the  qualitative  dye  mo¬ 
tion  observations  of  the  flow  patterns  could  not  be  confirm¬ 
ed  by  independent  pitot  tube  measurements.  Also  the  more 
modern  technique  of  laser  or  hot-wire  anemometry  have  not 
been  applied,  mainly  due  to  time  limitations. 

The  flow  visualisation  technique  used  has  introduced 
further  limitations,  in  that  secondary  flow  patterns  could 
only  be  observed  in  the  range  of  Reynolds  number  from  400 
to  1400  (Dean  numbers  676  to  2366) .  Below  this  range  the 
observations  were  made  impossible  by  the  dye  dispersion  and 
the  long  time  required  for  the  flow  stabilization,  and  above 
the  range  by  the  disturbances  produced  by  separation  zones. 

The  dye  method  also  required  the  continuous  use  of 
city  water, which  was  associated  with  presence  of  air 
bubbles,  visible  on  some  figures. 
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Fippre  1  Apparatus  Layout 
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Table  1 


Apparatus  Data 


Internal  diameter  of  the  pipe  and 

the  elbows 

d 

31  .75  mm 

Radius  of  curvature  of  pipe  centre¬ 
line 

R 

44.45  mm 

Dimensionless  curvature 

a/R 

0.357 

Relation  between  the  Reynolds  and 

Dean  numbers 

D  =  1.69  Re 

Spacer  lengths 

L 

0;  27  mm;  57  mm; 

100  mm 

Dimensionless  spacer  lengths 

L* 

0;  0.85;  1.8;  3- 

Inlet  pipe  length 

51  mm 

Outlet  pipe  length 

508  mm 

Inlet  tank  dimensions  (length  x 
width  x  height) 

305  mm  x  305  mm 

124  mm 

Outlet  tank  dimensions  (length  x 
width  x  depth) 

610  mm  x  202  mm 

203  mm 

Length  of  the  long  inlet  pipe 
(for  fully  developed  flow) 

b 

3664  mm 

Internal  diameter  of  the  pipe  as 

above 

A 

32.94  mm 

Effective  adapter  length 

28 . 6  mm 

Distance  from  the  water  level  in  the 

outlet  tank  to  the  outlet  pipe  axis 

approx.  75  mm 

Temperature  of  water 

22°  C 

Table  1 


continued 


r 

Kinematic  viscosity  of  water 

V 

0.95x10 6  m2/s 

Material  of  the  elbows  and  the  tanks 

plexiglass 

Diameter  of  the  pipes  injecting  the 

dye  mixture 

1  mm 

Camera  /  lens 

Canon  F  l/FD  55 

Close-up  lenses 

58  mm:  450;  240 

Film 

Kodak  Tri  -  X 

Lighting 

Fluorescent 
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Figure  3  Apparatus  -  General  View  (  Y  =  0°  ;  L'=  3.I5  ) 
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Figure  4  Apparatus  -  Elbows  (  ^  =  90°  ;  L*=  0  ) 


Figure  5  Apparatus  -  General  View  (  ¥  =  180°  ;  L*=  3*15  ) 
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Table  2 

Cross-Sections  and  Dye  Pipes  Positions 


X- 

s. 

DYE  PIPES* 

EXACT  POSITION 

NOTES 

A 

upstream  tank 

movable  -  any  position 

at  x-section  available 

Bl&Br 

1-st  elbow  inlet 

Cl&Cr 

3  mm  before 

1-st  elbow  outlet 

1 

1-st  elbow  outlet 

when  L'=  0,  X1=X2=X3 

i 

P 

_ 

spacer  ** 

rotated  with  spacer 

2 

spacer  mid-point 

when  L'=  0,  X1=X2=X3 

3 

2-nd  elbow  inlet 

Ki,Ko,Kl ,Kr 

13  mm  before 

2-nd  elbow  inlet 

4 

Ll,Lr 

2-nd  elbow  mid-point 

5 

2-nd  elbow  outlet 

Mi, Ml, Mr 

11  mm  behind 

2-nd  elbow  outlet 

6 

Nt 

42  mm  behind 

2-nd  elbow  outlet 

can  be  rotated  with  out- 

let  pipe  around  this  pi¬ 
pe  axis 

7 

OcC 

166  mm  behind 

2-nd  elbow  outlet 

*  for  the  dye  pipes  nomenclature  see  Figure  9 
**Dye  Pipes  in  Spacers  (distance  from  XI  in  mm): 


1 

DYE  PIPS 

r~ 

D 

E 

F 

G 

H 

I 

J 

3.15 

14 

25 

40 

52 

63 

7  8 

89 

SPACER  lmQ 

13 

21 

29 

38 

48 

LENGHT  0.85 

11 

19 

90*<  V<  180°  0*4V4  90' 
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Figure  6  Mirrors  and  Camera  Position 


CHAPTER  IV 


EXPERIMENT 

4.1  Experimental  Range 

As  can  be  concluded  from  Section  1.2,  no  references  on 
experimenbal  flow  pabbern  invesbigabions  in  bwo  successive 
pipe  elbows  or  bends  wibh  Reynolds  numbers  lower  bhan  10^ 
have  been  found.  In  bhe  presenb  work,  due  bo  bhe  apparabus 
limibabions  menbioned  in  Secbion  3.3,  only  bhe  lower  parb 
of  bhis  unreporbed  range  was  invesbigabed :  bhe  flow  pabbern 
invesbigabion  covered  bhe  range  of  bhe  Reynolds  numbers 
from  400  bo  1400  (Dean  numbers  from  646  bo  23 66)  and  bhe 
separabion  invesbigabion  covered  bhe  range  Re  from  400  bo 
2500  (D  from  676  bo  4225)  •  The  upper  limib  of  bhe  former 
was  resbricbed  by  bhe  burbulenb  disburbances ,  bub  bhe  upper 
limib  for  bhe  labber  was  exbended  bo  bhe  higher  flow  vel- 
ocibies  because  for  bhis  case  bhe  burbulence  did  nob  disburb 
bhe  observabions .  Selecbion  of  bhe  flow  velocibies  for  bhe 
illusbrabed  phobos  was  made  based  on  bhe  facbors  leading  bo 
bhe  good  qualiby  picbures,  bub  due  bo  bhe  weak  dependance 
of  bhe  flow  pabbems  on  bhe  Reynolds  (Dean)  number  over  bhe 
range  invesbigabed,  bhis  approach  seems  bo  be  correcb. 

To  provide  a  sabisfacbory  number  of  measuremenb  poinbs, 
bhe  posibion  of  bhe  separabion  lines  and  bhe  reabbachmenb 
lines  were  measured  for  8  or  9  differenb  flow  rabes. 

The  upper  limib  of  bhe  nexb  variable,  bhe  dimensionless 
spacer  lengbh  L' ,  was  chosen  equal  3*15  on  bhe  basis  of  ob¬ 
servabions  showing  no  observable  relabionship  bebween  bhe 
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flow  characteristics  and  the  spacer  length  for  values  of 
L*  larger  than  3.15*  For  the  purpose  of  the  separation 
position  measurements,  in  the  case  of  the  uniform  velocity 
profile , all  four  spacer  lenths  were  used.  The  experiments 
with  the  fully  developed  flow  at  the  inlet  required  fewer 
spacer  lengths.  The  measurements  using  only  two  of  them, 

L'  =  0  and  L'  =  3. 15, have  provided  enough  data  for  conclu¬ 
sions  and  comparisons. 

The  last  dimensionless  variable,  the  angle  was  va¬ 
ried  between  0°  and  180°  .  The  quantitative  results  of  the 
separation  measurements  were  obtained  for  the  values  of  the 
angle  V:  0*  ,  45°  ,  90°  *  135°  and  180®  for  the  constant  pres¬ 
sure  tank  and  for  15° ,  90°  and  180°  for  the  fully  developed 
flow  at  the  inlet. 

4.2  Experimental  Procedure 

The  procedure  of  the  experi¬ 
ment  was  separated  into  two  major  parts.  For  both  of  them 
the  water  temperature  was  precisely  adjusted  to  match  the 
surrounding  air  temperature  (22® C),  the  apparatus  was  sys¬ 
tematically  flushed  with  a  strong  jet  of  water  to  remove 
air  bubbles  and  the  flow  rate  was  coarsely  controlled  by  the 
valves  and  the  flowrator  meter,  and  precisely  measured  using 
the  flow  rate  measuring  tank  and  stopwatch.  Also,  in  the 
both  cases,  a  few  minutes  of  time  was  required  for  the  flow 
to  stabilize  after  any  change  of  the  flow  setting  (up  to 
10  min  for  the  lowest  flows) . 


The  first  procedure  was  applied  to  all  qualitative  ob¬ 
servations,  such  as  the  flow  pattern  investigation.  Con¬ 
ducting  this  procedure,  all  variable  parameters  were  changed 
within  the  ranges  described  in  the  previous  section  and  the 
dye  mixture  was  injected  at  different  dye  points  using  the 
trial -and-err or  method,  until  a  detailed  understanding  of 
the  observed  effects  had  been  reached.  Next,  the  flow  was 
photographed  for  certain  apparatus  settings,  chosen  to  ob¬ 
tain  the  best  photographic  results. 

A  different  procedure  was  applied  to  the  quantitative 
observations  of  the  separation  positions.  In  this  case  for 
the  apparatus  settings  given  in  the  previous  section,  the 
dye  mixture  was  injected  by  the  dye  pipe  located  in  the  area 
where  separation  probably  existed.  After  a  few  minutes,  re¬ 
quired  for  the  dye  stabilization,  the  distances  between  the 
extremal  upstream  dye  positions  and  the  reference  cross- 
sections,  and  L^,  were  measured  as  the  separation  lines 
positions.  The  reattachment  point  position  was  found  by 
looking  for  the  place  where  the  injected  dye  mixture  moves 
in  both  directions  simultaneously:  upstream  and  downstream. 

In  this  case  the  dye  pipes  located  in  a  spacer  downstream 
from  the  separation  line  position  were  used.  Due  to  the 
vortices  existing  behind  the  separation  zone,  especially  for 
the  higher  flow  velocities,  and  the  experimental  procedure 
employed,  the  numerical  results  obtained  for  the  reattach¬ 
ment  point  position  should  be  treated  as  approximated. 

All  numerical  results  were  recalculated  in  order  to 
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obtain  dimensionless  values,  L£,  L^, 


and  Ly 


4.3  Presentation  of  Results 

In  this  thesis  the  results  abe  presented  mainly  in  the 
form  of  graphs  and  tables  (numerical  results),  and  sketches 
and  photographs  (qualitative  results) .  Most  of  them  follow 
immediately  the  proper  section  of  the  text,  except  the 
tables  of  the  separation  position  measurements  which  are 
located  in  Appendix  for  reference. 

The  basis  for  the  interpretation  of  the  sketches  and 
photographs  are  Figures  7  &  8,  where  the  information  about 
the  cross-sections'  observation  points  and  the  cross-sections' 
locations  is  given  (for  the  exact  cross-sections'  locations 
see  also  Table  2) .  The  illustrative  material  also  contains 
the  data  of  the  dye  pipes  in  use  and  the  symbols  of 
the  flow  perpendicular  to  the  observer.  For  interpretation, 
see  Table  2  and  Figure  9*  Additionally,  the  following  notes 
should  be  kept  in  mind: 

-  the  velocities  v  and  w  are  at  least  one  order  of  mag¬ 
nitude  lower  than  the  average  flow  velocity  U 

-  in  certain  cases  when  the  dye  mixture  was  injected 
with  large  initial  momentum  to  investigate  the  flow  inside 
the  pipe,  the  jet  effect  must  be  taken  into  account  when 
the  photographs  are  studied 

-  in  practice,  it  was  impossible  to  obtain  pictures  of  the 
second  elbow  and  the  outlet  pipe  cross-sections  photographed 
from  the  cross-sections  observation  point  shown  in  Figure  8; 


' 

instead  the  camera  was  located  as  shown  in  Figure  6  or  be¬ 
hind  the  outlet  tank  (final  vortex  pictures)  which  means  the 
photographs  obtained  are  the  reverse,  around  the  vertical 
axis,  of  what  would  have  been  obtained  using  the  observa¬ 
tion  point  of  Figure  8 

-  the  "outer  curvature"  and  the  "inner  curvature"  of  the 
elbows  are  defined  in  Figure  8  (see  picture  in  the  upper 
right  corner)  . 
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Figure  7  Coordinate  System 
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FIG  8  NOTATION  OF  ELBOWS  AND  CROSS-SECTIONS 
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FIRST  ELBOW  SECOND  ELBOW 

DP's-B  and  C  DP's-K,L  and  M 


I .  C  .  1 


DPr 


DPo 


SPACERS  AND  OUTLET  PIPE 
DP's-Dot  to  J*  and  N*,  0^ 

0.C.2 


INLET  TANK,  POSITION 
SHOWN  GRAPHICALLY  (  DP -A) 


Note:  For  the  observation  positions  of  the  above 
cross-sections  see  Figure  8 


observation  point 


I 
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flow  direction 

®  perpendicular  flow  symbol 

'/////,  low  shear  stress  and  /  or  separation  zone 


Figure  9  Standard  Cross-Sections  Presentation,  Dye  Pipes 
Locations  at  Cross-Sections  and  Flow  Symbols 


CHAPTER  V 


SOME  OBSERVED  FLOW  PHENOMENA 
5.1  Secondary  Flow  and  Influence  .of  the  Second  Elbow 

The  general  character  of  fluid  motion  in  a  bend  is  well 
known  (Singh  197^) s  as  the  fluid  enters  the  bend  from  the 
reservoir,  a  boundary  layer  similar  to  that  in  a  straight 
pipe  develops  at  the  wall.  Thus  the  region  consists  of  two 
parts:  (i)  a  thin  layer  near  the  wall  where  the  viscous 
forces  are  balanced  by  the  inertia  forces  and  the  centrifu¬ 
gal  force  has  a  second-order  effect  near  the  inlet  and  (ii) 
an  inviscid  core  in  which  the  centrifugal  force  due  to  the 
circular  motion  of  the  main  body  of  the  fluid  along  the  pipe 
is  balanced  by  the  pressure  gradient  directed  towards  the 
centre  of  curvature.  However,  as  the  flow  develops  down¬ 
stream  the  bend,  the  curvature  effects  become  as  important 
as  the  viscous  effects  and  the  inertia  forces  in  the  boun¬ 
dary  layer,  which  in  turn  influence  the  flow  in  the  core 
due  to  fluid  particles  displacement.  The  secondary  flow  of 
the  fluid  due  to  curvature  consists  of  the  transverse  flow  of 
the  slower-moving  fluid  particles  in  the  boundary  layer 
from  the  outside  of  the  bend  towards  the  inside  and  the 
cross-flow  of  faster-moving  fluid  particles  from  inside  to 
the  outside  of  the  bend.  Thus,  the  secondary  boundary  layer 
behaves  as  a  reservior  receiving  the  fluid  moving  towards  the 
outer  curvature,  and  also  acts  as  a  source  because  of  the 
fluid  leaving  it  at  the  inner  curvature. 
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If  a  sufficiently  long  coiled  pipe  has  a  constant  radius 
and  a  constant  cross-section,  the  current  flow  eventual¬ 
ly  becomes  fully  developed,  that  is,  the  velocities  do  not 
vary  with  distance  along  the  pipe  axis.  Although,  in  view 
of  the  convective  mass  transfer  due  to  secondary  flow,  the 
developing  flow  requires  a  much  shorter  entrance  length  to 
become  fully  developed  in  comparison  to  a  straight  tube,  as 
was  found  by  Squire  ( 195*0  »  ‘the  transition  region  is  still 
relatively  long  and  the  flow  does  not  become  fully  develop¬ 
ed  for  large  dimensionless  curvature  90°  elbows. 

The  conclusion  obtained  from  the  above  consideration 
is  that  the  secondary  flow  starts  at  the  portion  of  the 
elbow  cross-section  where  the  curvature  is  the  largest 
(close  to  the  inner  curvature) .  When  two  successive  pipe 
elbows  are  located  to  form  "S"  shape  (V  =  180*  ) ,  the  influ¬ 
ence  of  the  second  elbow  on  the  flow  in  the  spacer  is 
strongest  close  to  its  inner  curvature,  making  the  stream¬ 
lines  in  the  spacer  almost  parallel  to  the  pipe  axis  on 
the  inner  side  closer  to  the  second  elbow  center  of  cur¬ 
vature  (Figure  10)  and  pushing  the  secondary  flow  circula¬ 
tion  to  the  outer  side. 

Another  conclusion  will  be  illustrated  in  Chapter  VI: 
when  the  flow  with  the  secondary  circulation  enters  a  second 
elbow,  the  transverse  flow  occurs  from  the  outside  of  the 
second  elbow  towards  the  inside.  This  flow  exists  in  a  thin 
layer  along  the  elbow  walls  and  is  independent  of  the  secon¬ 
dary  circulation  inside  the  flow  core. 
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5.2  Separation 

During  the  experiment,  the  phenomenon  of  boundary  layer 
separation  has  been  observed.  In  order  to  explain  this 
phenomenon,  let  us  consider  first  'the  flow  about  a  infinite¬ 
ly  long  circular  cylinder  (  Figure  11,  by  Schlichting,  1968  ). 
In  the  case  of  potential  flow,  for  particles  moving  along  DE, 
there  is  a  transformation  of  pressure  into  kinetic  energy, 
the  reverse  taking  place  along  EF,  so  that  a  particle  arri¬ 
ves  at  F  with  the  same  velocity  as  it  had  at  D.  In  the  actual 
case  of  a  viscous  flow,  the  boundary  layer  forms  and  a 
fluid  particle  which  moves  in  the  immediate  vicinity  of  the 
wall  in  the  boundary  layer  remains  under  the  influence  of 
the  same  pressure  field  as  that  existing  outside  and  impres¬ 
sed  on  the  b. 1.  However,  owing  to  the  large  friction  forces 
in  the  thin  boundary  layer  such  a  particle  consumes  so  much 
of  its  kinetic  energy  on  its  path  from  D  to  E  that  the  re¬ 
mainder  is  too  small  to  surmount  the  "pressure  hill"  from 
E  to  F.  Such  a  particle  can  not  move  far  into  the  region 
of  increasing  pressure  between  E  and  F  and  its  motion  is 
eventually  arrested.  The  external  pressure  causes  it  then 
to  move  in  the  opposite  direction.  This  reverse  motion 
gives  rise  to  a  vortex  whose  size  is  increased  until  it 
becomes  separated  and  moves  downstream  in  the  fluid.  Flow 
in  the  boundary  layer  near  a  point  of  separation  is  shown 
diagrammatically  in  Figure  11  -  B. 

The  above  example,  although  it  gives  the  explanation 
of  separation,  is  not  an  analogy  of  the  fluid  flow  inside 
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the  elbow.  To  make  an  approach  to  such  analogy,  con¬ 
sider  now  the  flow  in  the  duct  shown  in  Figure  12,  assuming 
it  is  infinitely  wide.  Such  a  duct  is  the  two  dimensional 
model  of  the  90°  large  dimensionless  curvature,  constant  ra¬ 
dius  elbow. 

Inside  this  elbow  model  the  kinetic  energy  of  a  particl 
at  point  E  is  relatively  high,  also  the  outer  wall  changes 
the  direction  of  the  flow  and  keeps  the  point  of  maximum 
velocity  closer  to  the  center  of  curvature.  In  this  situa¬ 
tion  the  separation  occurs  due  to  adverse  pressure  gradient 
existing  between  the  inner  curvature  and  the  area  downstream 
from  the  elbow .  Inside  the  elbow  the  centrifugal  force 
creates  a  pressure  difference  between  the  inner  and  the  out- 
ter  curvature  of  the  elbow,  with  the  lower  pressure  closer 
to  the  inner  curvature.  This  moves  the  maximum  velocity 
point  outward  and  the  flow  at  the  inner  curvature  decelera¬ 
tes  .  Downstream  from  the  elbow,  the  point  of  the  maximum 
velocity  moves  again  to  the  center  of  the  pipe  and  pressure 
measured  at  the  inner  curvature,  along  the  pipe,  increases. 
The  separation  ends  with  the  reattachment  point  when  pres¬ 
sure  starts  to  decrease  again. 

The  two-dimensional  elbows  can  be  combined  in  the  two 
basic  configurations:  "U"  shape  and  "S"  shape  with  or  with¬ 
out  a  spacer  between  them.  The  "S"  configuration,  especial¬ 
ly  without  or  with  a  short  spacer,  gives  an  interesting 
contribution  to  the  above  considerations:  the  centrifugal 
force  inside  the  second  elbow  causes  a  relatively  high  pres- 


■ 

sure  increase  in  the  flow  between  the  inner  curvature  of 
the  first  elbow  and  the  outer  curvature  of  the  second  one. 
This  pressure  gradient  promotes  the  separation  and  as  ex¬ 
pected  moves  the  separation  line  upstream.  The  strength  of 
the  centrifugal  force  depends  on  the  flow  velocity,  there¬ 
fore  the  length  of  the  separation  zone  also  depends  on  it. 

On  the  other  hand,  for  the  "U"  configuration  without  spacer 
the  separation  is  expected  to  be  reduced  due  to  the  direc¬ 
tion  of  the  second  elbow  centrifugal  force.  This  is  re¬ 
lated  t o  the  "first"  separation,  which  originates  in  the 
first  elbow,  or  the  spacer,  or  the  second  elbow  when  there  i 
no  spacer,  but  a  separation  also  can  be  expected  in  the  se¬ 
cond  elbow  and  the  outlet  pipe, called  the  second  separation. 

The  described  phenomena  of  the  first  and  second  sepa¬ 
ration  and  the  influence  of  the  flow  velocity,  are  even 
better  expressed  for  the  three-dimensional  case  of  the  el¬ 
bows  with  the  circular  cross-section  (Figures  1 3 »  and 

15) t  due  to  the  secondary  circulation  which  is  changing 
the  character  of  the  flow.  Instead  of  the  steady  circula¬ 
tory  motion  inside  the  separation  zones,  as  in  the  two-dim¬ 
ensional  case,  the  presence  of  the  secondary  flow  can  create 
backflow  supplied  by  the  secondary  circulation  boundary  lay¬ 
er  flow  and  discharged  in  the  plane  dividing  two  vortices 
of  the  secondary  flow,  as  shown  schematically  in  Figure  16 
(see  also  Figures  17  through  20  inclusive).  This  backflow 
has  been  observed  over  the  whole  investigated  range  of  the 
flow  in  the  first  elbow,  where  the  two-coil  structure  of  the 
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secondary  flow  is  always  well-defined.  The  phenomenon  also 
occurs  in  the  second  elbow  for  Re>l400,  where  the  second¬ 
ary  flow  structure  depends  on  a  mutual  position  of  the  el¬ 
bows  (Figure  21 ) . 

The  results  of  the  separation  measurements  are  graph¬ 
ically  presented  in  Figures  22  through  39  inclusive.  For 
the  constant  pressure  tank  the  first  separation  line  posi¬ 
tion  is  shown  in  Figures  22  through  25,  the  first  separa¬ 
tion  reattachement  line  position  in  Figures  26  through  29 
and  the  second  separation  line  position  in  Figures  30  through 
33*  The  same  measurements  for  the  fully  developed  flow  for 
selected  apparatus  settings  are  presented  in  Figures  3^ 
through  39* 

These  graphs  (for  numerical  results  see  Appendix)  show 
the  similar  relationships  as  previously  concluded  to  the 
two-dimensional  case  observations.  For  angle  V  =  180°  the 
high  pressure  region  present  on  the  outer  side  of  the  flow 
in  the  second  elbow  moves  the  first  separation  line  up¬ 
stream.  When  angle  V<  180®  this  effect  is  weaker;  when 
V  =  0*  ,  the  low  pressure  region  present  on  the  inner  side 

of  the  flow  in  the  second  elbow  moves  the  separation  line 
downstream  and  for  the  case  without  spacer  (L'  =*  0)  prevents 
the  separation  from  occurring  in  the  first  elbow. 

The  influence  of  the  angle  V  on  the  first  separation 
line  position  is  observed  even  for  the  longest  spacer 
(L'  =  3.15),  but  the  influence  of  the  spacer  length  is  weak 
for  L'>  1.8  due  to  the  weaker  influence  of  the  second  elbow 


on  the  flow  upstream. 

The  first  separation  reattachment  line  position  is 
weakly  affected  "by  the  angle  V  -  the  curves  shown  in  Figure 
26  through  29  inclusive  are  located  close  one  to  another 
except  the  curve  plotted  for  V  =  180*  in  Figure  26.  This 
one  is  located  more  downstream,  due  to  the  action  of  the 
backflow,  already  described.  The  influence  of  the  spacer 
length  is  observed  only  for  L'  =  0  and  L'  =  0.85,  where  the 
presence  of  the  second  elbow  limits  the  first  separation. 

In  other  cases,  L' >  0.85,  the  length  of  the  separation, 
measured  between  the  cross-section  number  1  and  the  reat¬ 
tachment  line  ,  is  almost  the  same . 

The  second  separation  line  position  as  a  function  of 
the  angle  V  and  the  spacer  length  L'  can  be  probably  ex¬ 
plained  in  terms  of  the  secondary  flow  patterns.  In  the 
case  of  the  "U"  configuration  (V  =  0*  ) ,  the  pressures 
measured  perpendicularly  to  the  axial  velocity  direction 
between  a  point  located  at  the  inner  curvature  of  the  sec¬ 
ond  elbow  and  a  point  located  directly  downstream  on  the 
wall  of  the  outlet  pipe,  should  show  relatively  small  dif¬ 
ference  with  the  lower  pressure  being  in  the  elbow.  This 
would  be  due  to  the  relatively  strong  secondary  cross-flow 
at  the  outlet  of  the  second  elbow,  directed  from  the  inside 
to  the  outside  of  the  elbow  (see  Section  5*1) •  Finally  the 
adverse  pressure  gradient  is  weak  and  the  separation  is  re¬ 
latively  small.  In  the  case  of  the  "S"  conf iguration 
(V  =  180°  )  the  secondary  circulation  in  the  outlet  pipe  is 
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more  complicated  and  forms  a  four  coil  structure  featuring 
the  two  upper  coils  created  by  the  first  elbow,  and  the  two 
lower  coils  created  by  the  second  one.  The  pressures 
measured  perpendicularly  to  the  flow  direction  between  a 
point  located  at  the  inner  curvature  of  the  second  elbow 
and  a  point  located  directly  downstream  on  the  wall  of  the 
outlet  pipe,  should  show  significant  difference  with  the 
lower  pressure  being  in  the  elbow.  The  adverse  pressure 
gradient  is  strong  and  the  separation  line  is  located  up¬ 
stream,  relative  to  the  case  of  "U"  configuration.  For  any 
other  values  of  the  angle  the  separation  line  occupies 
intermediate  positions. 

The  relationship  between  the  second  separation  line 
position  and  the  spacer  length  L'  is  probably  associated 
with  the  strength  of  the  two  upper  vortices  created  by  the 
first  elbow.  These  vortices  in  the  second  elbow  are  weaker, 
when  the  spacer  is  longer,  which  allows  the  second  elbow  to 
develop  stronger  lower  vortices  and  reduces  the  adverse 
pressure  gradient. 

The  influence  of  the  Reynolds  (Dean)  number  has  been 
observed  as  the  rapid  change  of  the  separation  lines  posi¬ 
tions,  which  is  associated  with  the  transition  from  laminar 
to  turbulent  flow  in  the  boundary  layer.  This  position 
change  is  well-defined  for  the  first  separation  and  weaker 
for  the  second  one.  The  possible  explanation  is  that  the 
flow  approaching  the  second  separation  zone  is  already  par¬ 
tially  turbulent.  The  relationship  between  the  reattachment 
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point  position  and.  the  Reynolds  (Dean)  number  also  shows, 
for  the  spacers  lengths  L*  =  0.85  and  L'  =  1.8,  a  similar 
rapid  change  of  the  position. 

The  fully  developed  velocity  profile  at  the  elbows  en¬ 
trance  has  introduced  significant  differences  (compared  with 
the  uniform  velocity  profile)  for  the  first  separation  line 
position  only.  Due  to  the  lower  velocity  near  the  wall, 
this  type  of  flow  is  more  sensitive  to  the  pressure  changes 
and,  for  example,  for  the  "S"  configuration  the  first  se¬ 
paration  line  position  moves  upstream  very  rapidly  as  the 
Reynolds  number  is  increased. 

5.3  Turbulence 

The  fact  that  separation  in  steady  flow  occurs  only 
in  decelerated  flow,  leads  directly  to  the  conclusion  that 
there  must  exist  a  point  of  inflection  of  the  velocity  pro¬ 
file  in  the  boundary  layer  (Schlichting  1968)  .  The  exis¬ 
tence  of  such  a  point  is  important  for  the  transition  from 
laminar  to  turbulent  flow.  It  provides  a  necessary  and 
sufficient  condition  for  the  amplification  of  disturbances, 
as  has  been  deduced  from  the  frictionless  Orr-Sommerf eld 
(stability)  equation.  This  equation  was  obtained  by  omit¬ 
ting  the  viscous  terms  in  the  full  Orr-Sommerf eld  equation, 
as  compared  with  the  inertia  terms, and  the  influence  of  vis¬ 
cosity  changes  the  preceding  conclusion  only  very  slightly. 
The  inflection  point  stability  criterion  means  that  a  fa¬ 
vourable  pressure  gradient  stabilizes  the  flow,  whereas  an 
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adverse  pressure  gradient  enhances  instability;  this  state¬ 
ment  helps  to  explain  the  basic  difference  between  White's 
(1929)  experiment,  where  the  critical  value  of  the  Reynolds 
number  for  transition  to  turbulence  was  found  to  be  Re=6000, 
and  the  present  experiment,  where  the  first  turbulent  in¬ 
stabilities  were  observed  for  Re=l407»  The  White  experiment 
v/as  conducted  for  a  long  pipe  of  relatively  small  dimension¬ 
less  curvature,  a/R=0.02,  where  no  separation  was  observed 
and  the  centrifugal  force  efficiently  stabilizes  the  flow. 

In  the  present  experiment  the  flow  was  observed  to  be  fully 
laminar  for  Re=796  (Figure  40)  and  for  Re=938  (Figure  41) 
when  the  laminar  vortices  produced  by  the  first  separation 
were  visible.  For  Re=l407  the  flow  center  remained  fully 
laminar  (figure  42)  but  the  first  separation  started  to  pro¬ 
duce  instabilities  which  are  probably  turbulent  (Figure  43) . 
This  is  based  on  the  data  from  Figure  22  through  25  inclu¬ 
sive,  showing  a  rapid  change  of  the  first  separation  line 
position  which  is  characteristic  for  the  boundary  layer 
transition  from  laminar  to  turbulent  conditions.  This  boun¬ 
dary  layer  transition,  when  the  center  of  the  flow  still  re¬ 
mains  laminar,  is  likely  to  be  associated  with  the  presence 
of  the  separation  in  the  entrance  cross-section  of  the  inlet 
pipe.  For  the  higher  Reynolds  numbers  the  laminar  vortices 
also  occur  at  the  outer  curvature  of  the  first  elbow 
(Figure  44)  and  instabilities  are  produced  by  the  second 
separation  (Figure  45) .  Further  increase  of  the  Reynolds 
number  shows  that  the  flow  center  still  remains  laminar 
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(Figure  46),  but  outside  the  laminar  center  the  instabili¬ 
ties  are  growing  (Fi^ire  47)  . 

Although,  as  has  been  shown,  White's  experiment  may 
not  be  directly  compared  with  the .present  ones,  both  of  them 
have  given  almost  the  same  Dean  numbers  for  the  transition 
to  turbulence,  if  the  instabilities  observed  in  this  ex¬ 
periment  for  Re=l407  are  assumed  to  be  turbulent.  These 
Dean  numbers  are  2400  and  2378,  respectively. 

The  results  obtained  also  show  no  contradiction  with 
the  iJerem  et  al  (1972)  measurements,  showing  only  the  la¬ 
minar  j.low  up  to  Re=5000  in  the  aortic  arch.  In  the  present 
experiment,  the  presence  of  the  laminar  center  for  Re=2221 
confirms  that  result,  but  further  comparisons  are  impossible 
due  to  the  different  entrance  conditions.  The  sharp  edge 
of  the  inlet  pipe  in  the  inlet  tank  (see  Figures  4  and  6) 
may  not  be  used  as  a  model  of  the  inlet  to  the  aortic  arch, 
due  to  the  ring-shaped  separation,  which  appears  inside  the 
inlet  pipe,  immediately  downstream  from  the  edge.  This 
separation  produces  observable  instabilities  which  disturbs 
the  flow  downstream.  This  is  so  even  for  "U"  configuration, 
when  there  is  no  first  separation  (Figures  48  and  49) . 


. 
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Figure  10 


Influence  of  the  Second  Elbow  on  the  Flow  in  the 
Spacer 

(  left  fig.:  Re  =  767,  V  =  180',  L’=  3.I5,  DP-AO) 
(  right  fig.jRe  =  767,  V  =  180°,  L’  =  3-15,  DP-EO  ) 
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•Figure  11  Separation  and  Vortex  Formation  on  a  Circular 
Cylinder  and  Flow  in  the  Boundary  Layer  Near 
a  Point  of  Separation  (S  -  point  of  separation) 


.Figure  12  Separation  in  a  Two  -  Dimensional  Elbow 
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Figure  13  First  Separation  as  a  Function  of  Reynolds  Number 
(  upper  fig.:  Re  =  400,  V  =  180*  ,  L'=  0,  DP  -  Cl  ) 
(  lower  fig.:  Re  =  650,  V  =  180°  ,  L'=  0,  DP  -  Cl  ) 
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Figure  14  Vortices  Due  to  the  Second  Separation 

(  Re  =  1055,  V  =  180® ,  L' =  0,  DP-A  and  Mi, 
aT  =  0.33s  ) 
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Figure  15  Second  Separation  -  Dye  Motion 

(  Re  =  403,  V  =  180°,  L’=  0,  DP -Mi,  a  T  =  2s  ) 
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Figure  16  Backflow  Streamlines  in  the  First  Separation 
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Figure  1?  Backflow  in  the  Spacer 

(  Re  =  938,  V  =  180°,  L*  =  3-15*  upper  fig.  DP -DO*, 
lower  fig.  DP-E200 ) 
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Figure  18 


First  Separation  -  Dye  Motion  Due  to  the  Backflow 
(  see  mirror  image  on  the  right,  Re  =  403>  =180° 

L '  =  0,  DP-A©  ,  AT  =  0.5s  ) 
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Figure  19  Backflow  in  the  First  Elbow 

(  Re  =  420,  W  =  180°  ,  L*=  0,  DP-Ko  ) 


Figure  20  First  Separation  -  Reattachment  Line 

(  Re  =  414,  V  =  180°,  L'  =  3-15,  DP-E0° ) 
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Figure  21  Second  Separation  as  a  Function  of  Reynolds 
Number  (  V  =  0°  ,  L’=  0,  DP-Mi,  upper  fig. 

Re  =  1530,  lower  fig. Re  =  2400  ) 


not  available 
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igure  22  First  Separation  Line  Position  versus  Reynolds  (Dean)  Number  as 
a  Function  of  Angle  t  for  a  Uniform  Velocity  Profile  and  L*  =  0 
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Figure  26  First  Separation  -  Reattachment  Line  Position  versus  Reynolds  (Dean) 
Number  as  a  Function  of  Angle  *  for  a  Uniform  Velocity  Profile  and  I/  =  0 
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Figure  <^8  First  Separation  -  Reattachment  Line  Position  versus  Reynolds  (Dean) 
Number  as  a  Function  of  Angle  T  for  a  Uniform  Velocity  Profile  and  L*  =  1.8 
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Figure  32  Second  Separation  Line  Position  versus  Reynolds  (Dean)  Number  as 
a  Function  of  Angle  V  for  a  Uniform  Velocity  Profile  and  L'  =  1. 


67 


15° not  observed 


68 


o  o 
o  o 
tH  o 


GJ 

pc; 


U-\ 

CM 


o 

o  O 

o  cc 

On  »-h 


II  II  II 

>  3-  *- 
•  >  a 


--3- 


o 

CM 


- 


Figure  3^  First  Separation  Line  Position  versus  Reynolds  (Dean)  Number  as 
Function  of  Angle  V  for  a  Fully  Developed  Velocity  Profile  and  L*  =  0 
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Figure  36  First  Separation  -  Reattachment  Line  Position  versus  Reynolds  (Dean)  Num¬ 
ber  as  a  Function  of  Angle  V  for  a  Fully  Developed  Velocity  Profile  and  L'  =  0 
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Figure  37  First  Separation  -  Reattachment  Line  Position  versus  Reynolds  (Dean) 

Number  as  a  Function  of  Angle  V  for  a  Fully  Developed  Velocity  Profile  and  L'  =3.15 
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igure  39  Second  Separation  Line  Position  versus  Reynolds  (Dean)  Number  as  a 
Function  of  Angle  tf  for  a  Fully  Developed  Velocity  Profile  and  L'  =  3.15 
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Figure  40  Fully  Laminar  Flow  with  the  Backflow 
(  Re  =  796,  V  =  180®  ,  L*  =  0,  DP-AO) 


Figure  41  Laminar  Vortces  Produced  by  the  First  Separation 
(  Re  =  938,  V  =  180°  ,  L '  =  0,  DP-AO) 
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Figure  42  Laminar  Flow  Center 

(  Re  =  1407,  V  =  180*  ,  L '  =  0,  DP-AO) 


Figure  43  Instabilities  Produced  by  the  First  Separation 
(  Re  =  1407,  V  =  180°  ,  L'=  0,  DP-A®) 
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Figure  44  Laminar  Disturbances  at  the  Outer  Curvature  of 

the  First  Elbow  (Re=1835,  V  =  180°,  L'=  0,DP-AQ) 


Figure  45  Instabilities  Produced  by  the  Second  Separation 
(  Re  =  1835,  *  =  180°,  L'  =  0,  DP-A©) 
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Figure  46  Laminar  Flow  Center 

(  Re  =  2221,  V  =  180°  ,  L'=  0,  DP-A  0 ) 


Figure  47  Instabilities  Outside  the  Laminar  Center 
(  Re  =  2221,  V  =  180°,  L‘=  0,  DP-A©) 


Figure  48  "U"  Configuration  Instabilities 

(  Re  =  1507,  V  =  0°  ,  L '  =  0,  DP-A®) 


Figure  49  "U"  Configuration  Instabilities 

(  Re  =  2200,  V  =  0°  ,  L'=  0,  DP-A®) 


CHAPTER  VI 


SECONDARY  FLOW  PATTERNS 
6.1  Introductory  Remarks 

The  purpose  of  "this  chapter  is  to  present  the  results 
of  the  qualitative  observations  of  the  secondary  flow  pat¬ 
terns.  Although  these  patterns  depend  strongly  on  the  spacer 
length  and  the  angle  V,  the  following  remarks  are  valid  for 
all  flow  types,  presented  in  Tables  3  and  4. 

nxcept  for  a  very  few  cases,  there  is  no  observable  in¬ 
fluence  of  the  Reynolds  (Dean)  number  on  the  secondary  flow 
patterns . 

-  The  secondary  flow  velocities  were  found  by  the  dye  motion 
observations  to  be  at  least  one  order  of  magnitude  smaller 
than  the  main  flow  velocity  U. 

-  The  secondary  flow  always  occurs  first  at  the  inner  side 
of  the  elbow  cross-section. 

-  The  secondary  flow  in  the  first  elbow  always  forms  the 
symmetrical  two  coil  structure,  as  described  in  Section  5*1* 
In  certain  cases  (low  Reynolds  number)  this  structure  may 

be  disturbed  by  the  influence  of  the  second  elbow,  but 
the  effect  is  not  significant. 

-  Regardless  of  the  flow  in  the  core,  in  the  thin  layer  a- 
long  the  walls  of  the  second  elbow  there  always  exists  flow 
directed  inward,  called  the  boundary  layer  flow  of  the  sec¬ 
ond  elbow  secondary  circulation.  The  velocity  of  this  flow, 
compared  with  the  mean  velocity  of  the  flow»U,  is  very  slow 
and  the  layer  is  very  thin,  therefore  even  a  minimal  outward 
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motion  inside  the  flow  center  is  sufficient  to  supply 
enough  fluid  to  meet  continuity  requirements.  This  flow 
has  been  observed  in  the  elbow  only. 

-  The  second  elbow  magnifies  the  secondary  flow  patterns 
generated  by  the  first  elbow  when  V  =  0°.  When  enhancement 
is  impossible,  for  example  when  the  first  elbow  secondary 
flow  patterns  don't  match  the  second  elbow  plane  of  cur¬ 
vature,  the  second  elbow  produces  its  own  secondary  flow 
patterns  starting  from  the  inner  curvature. 

-  In  the  outlet  pipe,  downstream  from  the  second  elbow,  for 
all  practical  purposes,  only  inertia  and  viscosity  forces 
are  acting  on  the  secondary  flow.  On  the  distance  not  lar¬ 
ger  than  approximately  10  diameters,  these  forces  change  the 
structures  of  the  two,  three  or  four  vortices,  present  in 
the  cross-section  number  5>  as  follows: 

*  the  smaller  vortices  combine  with  the  larger  one,  if 
they  are  in  the  same  direction, 

*  the  stronger  vortex  overwhelms  the  weaker  ones  if  are  in 
the  opposite  directions,  and  in  addition  the  weaker  ones 
are  more  quickly  damped  by  the  viscous  force, 

*  the  two  or  four  vortex  structures  remains  unchanged  in 
the  flow,  when  the  neighboring  vortices  are  in  opposite  di¬ 
rections  and  equally  strong  they  are  equally  damped  by  the 
viscous  force  and  have  been  observed  even  in  the  outlet  tank. 

-  Due  to  the  apparatus  construction,  results  have  been  ob¬ 
tained  for  the  values  of  angle  0°^  180®  ,  but  symmetry  al¬ 

lows  the  extrapolation  of  these  results  to  cover  180%  V  4  360° 
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6.2  "U"  Configuration  (¥  =  0°) 

This  configuration  (angle  V  =  0°)  produces  the  second¬ 
ary  flow  type  AO,  presented  in  Tables  3  and  4.  The  second 
elbow  amplifies  the  two  coil  structure  which  has  been  ob¬ 
served  as  far  as  the  outlet  tank.  Two  additional  outer- 
vortices,  reported  by  Rowe  (1970)  for  high  Reynolds  number 
flow,  Re  =  2.36  x  10-^  have  not  been  observed  (Figures  50, 

51  and  52 ) . 

6.3  "S"  Configuration  (V  =  180°) 

This  configuration  produces  secondary  flow  types  E0  or 
El  depending  upon  the  flow  velocity.  In  this  case,  the  sec¬ 
ondary  circulation  of  the  second  elbow  starts  at  the  inner 
curvature  of  the  second  elbow,  and  at  the  outlet  of  this 
elbow  the  final,  symmetrical  four  coil  structure  is  observed. 
This  structure  consists  of  two  upper  vortices  which  are  pro¬ 
duced  by  the  first  elbow  and  two  lower  ones  which  are  pro¬ 
duced  by  the  second  elbow,  and  persists  along  the  outlet 
pipe  when  the  Reynolds  number  is  higher  than,  approximately, 
450.  For  lower  values  of  Re,  the  vortices  produced  by  the 
second  elbow  are  weak  and  disappears  before  the  cross-sec¬ 
tion  number  7* 

The  presence  of  a  spacer  allows  the  second  elbow  to 
influence  the  flow  between  the  elbows.  This  influence  moves 
the  two  vortices  produced  by  the  first  elbow  to  the  spacer 
wall  located  directly  upstream  from  the  outer  curvature  of 
the  second  elbow  (see  points  "0"  in  Figure  8).  At  this 
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time,  the  streamlines  located  closer  to  the  opposite  side 
o!  a  spacer  and  upstream  from  the  inner  curvature  of 
ohe  second  elbow,  are  almost  parallel  (points  "I"  in  the 
spacer  in  Figure  8) .  This  area  without  vortices  grows, 
when  proceeding  from  the  cross-section  1  to  the  cross-sec¬ 
tion  3,  where  the  two  vortices  of  the  second  elbow  appear. 
Figures  54  through  57  inclusive. 

6.4  Flow  Patterns  for  0°<  180° 

The  three-dimensional  configuration  of  the  elbows 
centerline  allows  the  observation  of  more  complicated  flow 
patterns  than  described  in  Sections  6.2  and  6.3.  Slowly 
changing  the  angle  V  from  0*  to  180*  and  observing  the  sec¬ 
ondary  flow  behavior  between  the  elbows,  for  example,  in 
the  spacer  in  Figure  58,  it  can  be  seen  that  up  to  a  certain 
critical  value  of  the  angle  V,  the  flow  maintains  the  two 
coil  structure,  characteristic  for  the  "U"  configuration. 

To  do  this,  the  flow  acts  this  ways  the  plane  of  sym¬ 
metry  of  the  two  coil  structure  is  twisted  spirally  along 
the  flow  axis,  in  order  to  match  the  second  elbow  plane  of 
curvature,  where  the  centrifugal  force  is  acting.  This  type 
of  flow,  associated  with  the  "U"  configuration,  is  called 
"A"  (A1  -  A3  when  V  4  0^  and  is  discussed  later. 

The  opposite  situation  is  observed,  when  the  angle  V 
is  slowly  reduced  from  180*  to  0°  .  This  time,  down  to  a 

certain,  critical  value  of  the  angle  V,  the  flow  maintains 
the  four  coil  structure,  characteristic  for  the  "3" 


conf igu- 
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ration.  The  flow  behaves  such  that  the  plane  of  symmetry 
of  the  first  elbow  coils  (Figure  58)  is  spirally  twisted  to 
match  the  four  coil  structure  observed  at  the  outlet  of  the 
second  elbow  when  the  angle  V  =  180°  .  This  type  of  flow, 

associated  with  the  "S"  configuration  is  called  "B"  type 
(^2  when  V  *  180e  )  and  the  plane  of  symmetry  is  twisted  in 
opposite  direction  than  for  flow  type  "A" .  The  criti¬ 
cal  values  of  the  angle  t  for  the  flow  type  A  and  the  flow 
type  B  are  not  equal  and  f  crit  A<  *crit  A  when  0°^  * 4  180° 
or,  as  can  be  concluded  *  crit  A  >  *  crit  £  when  180°4*  4  3600  . 
Additionally,  there  is  no  direct  switch  from  the  type  A  to 
the  type  B,  except  in  the  case  when  L'=  0;  the  intermediate 
type,  called  "AB"  has  been  observed,  featuring  no  twist  of 
the  plane  of  the  two  coil  structure  symmetry  in  the  spacer, 
(Figures  58  and  from  69  to  71  inclusive).  This  type  has 
been  observed  to  be  unstable,  expecially  for  small  L'  values 
and  readily  changes  to  either  type  A  or  type  E.  The  area 
between  the  *  cr^  A  and  T  crq-t  r  is  ibe  same  where  the  flow 
type  AB  has  been  observed.  The  critical  values  of  the  angle 
^  are  likely  to  depend  strongly  on  random,  external  vibra¬ 
tions  and  the  rate  at  which  the  angle  V  is  changed.  There¬ 
fore,  results  of  the  measurements  have  had  a  significant  scat¬ 
ter  and  in  Figure  59  the  more  consistent  data  have  been 
collected . 

As  has  been  mentioned  previously  for  the  flow  patterns 
type  A  and  B,  the  flow  in  the  spacer  or  in  the  second  elbow 
when  there  is  no  spacer,  acts  such  that  the  plane,  div- 
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viding  the  two  symmetrical  vortices  produced  by  the  first 
elbow,  is  twisted  in  order  to  match  it  to  the  second  elbow 
plane  of  curvature.  Eut  such  twist,  expecially  when  L'  is 
small,  disturbs  the  symmetry  of  the  two  vortices,  making 
one  of  them  stronger  and  the  second-weaker,  adds  a  rotation 
to  the  whole  structure  and  the  twist  can  not  be  completed 
before  the  second  elbow  entrance.  Finally,  under  the  forces 
in  the  second  elbow,  the  flow  patterns  are  still  changing 
downstream  from  this  elbow  entrance  (see  Tables  3  and  4). 

An  interesting  confirmation  of  the  above  can  be  given 
by  the  observation  of  the  flow  type  A  when  there  is  no  spacer 
between  the  elbows.  For  such  case  the  flow  type  A  has 
been  observed  up  to  V  =  165°  .  In  the  range  of  the 
0C<  V  £  I350  the  two  coil  structure  of  the  first  elbow  is 
twisted  inside  the  second  elbow,  as  described  previously 
and  this  is  the  version  A1  (Tables  3  and  4,  Figures  60,  61 
and  62).  When  the  angle  t  is  larger,  approximately  120°  to 
lo0°  ,  some  of  the  fluid  particles  are  not  able  to  make  such 
a  long  turn,  due  to  the  presence  of  the  backflow  zone  on 
the  outer  side  of  the  second  elbow  (Section  5*2).  These 
particles  make  the  smaller  turn,  dividing  the  larger  coil 
in  two,  as  shown  in  Table  4  and  Figures  63,  65  and  67).  For 
the  largest  ¥  allowable  for  the  flow  type  A  (from  155°  "to 
I650)  the  two  vortices  obtained  from  the  larger  one,  as  descri¬ 
bed  above  ,  are  smaller.  The  additional  space  in  the  fifth 
cross-section  is  occupied  then  by  the  small  vortex  created 
by  the  shear  stress  (Table  4,  Figures  64,66  and  68).  The 
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<^^^'^'erences  between  the  two  versions  of  flow  (called  A2  and 
A3  respectively)  are  shown  in  Figures  65  through  68  inclusive 
The  versions  A2  and  A3  have  not  been  observed  for  spacer 
lengths  greater  than  L’  =  0.85  and  for  the  Reynolds  number 
Re<450  which  confirms  the  association  of  these  versions  of 
flow  with  the  backflow  in  the  second  elbow. 

xhe  details  of  the  flows  types  AB  and  B2  are  shown  in 
Figures  69  through  77  inclusive. 

One  of  the  most  important  features  of  the  flows  des- 
^ ^ ® ^  this  section  is  the  final  effect  observed  down- 
otream  from  the  cross-section  number  7 •  This  effect  is  al¬ 
ways  (qualitatively)  the  same,  regardless  of  the  angle  'f  and 
the  type  of  the  flow  inside  the  second  elbow,  and  has  the 
form  of  the  single  vortex,  turning  counterclockwise,  when 
0  <  t  <  180  and,  concluding  by  the  symmetry,  clockwise,  when 
180  <360  .  This  is  shorn  in  Figure  78  (due  to  the  obser¬ 

vation  point,  the  motion  on  the  photographs  is  clockwise), 
where  all  pictures  show  the  same  direction  of  rotation, 
except  one  of  the  pictures  taken  for  f  =  30°  ,  where  the 
second  vortex  (the  flow  is  type  Al)  has  been  too  strong  to 
be  overwhelmed  or  damped  by  the  viscosity. 
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Table  4  Secondary  Flow  Patterns 

Flow  Types  with  a  Short  Spacer  (  L'^  0.85  ) 
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Figure  50  Flow  Type  AO  -  Left  Coil 

(  Re  =  403,  V  =  0°  ,  L'=  0,  DP-Ko  ) 
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Figure  51  Flow  Type  AO  -  Two  Coil  Structure 
(  Re  =  403,  V  =  0°  ,  L'=  0,  DP-AO) 


1-st  elbow 
inlet 


1-st  elbow 
mid-poin* 


Figure  52  Flow  Type  AO  -  Left  Coil 

(  Re  =  403, y  =  0°,  L’=  0,  DP-AO) 
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Figure  53  Flow  Type  BO  -  Lower  Coils 

(  Re  =  681,  V  =  180*  ,  L'  =  3-15,  DP -HO0-) 
low  initial  momentum) 
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Figure  5^  Flow  Type  BO  -  Upper  Coils 

(  Re  =  812,  V  =  180°  ,  L'=  3-15.  DP-Ko-high  ini¬ 
tial  momentum  ) 
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Figure  55  Flow  Type  B1  -  Upper  Coils 

(  Re  =  380,  V  =  180°,  L '  =  0,  DP-Ko-high 
initial  momentum  ) 
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Figure  56  Flow  Type  Bl  -  Lower  Coils 

(Re  =  380,  V=  180*  ,  L ' =  3.15,  DP-I0°) 
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Figure  57  Flow  Type  B1  -  Right  Coils 

(  Re  =  380,  ^  =  180°,  L‘=  3. 15,  DP-I0° ) 
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Figure  58  Basic  Flow  Types 

(  Re  =  422,  L'  =  3. 15,  DP-A  © ,  left  fig.  V  =  90°  , 
central  fig.  V  =  650 ,  right  fig.  V  =  90°  ) 
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Figure  59  Area  of  the  Flow  Types  Existence 
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Figure  60  Flow  Type  A1  -  Left  Coil 

(  Re  =  659,  V  =  135* ,  L'=  0,  DP-C1  ) 
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Figure  6l  Flow  Type  A1  -  Two  Coil  Structure 

(  Re  =  485,  V  =  45°,  L'=  3. 15,  DP-A  O  ) 
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Figure  62  Flow  Type  A1  -  Left  Coil  and  Final 
(  Re  =  659,  V  =  90°  ,  L'  =  0,  DP-Cr  ) 


Vortex 
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Figure  63  Flow  Type  A2  -  Lower  Left  Coil  and  Backflow 
(  Re  =  659,  V  =  160° ,  L'=  0,  DP -Cl  ) 


Figure  64  Flow  Type  A3  -  Upper  Right  and  Lower  Left  Coils 
(  Re  =  508,  V  =  150°  ,  L ' =  0,  DP-Cr  ) 
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Figure  65  Flow  Type  A2  -  Upper  and  Lower  Right  Coils,  Final 
Vortex  (Re  =  508,  f  =  120* ,  L’=  0,  DP-Lr-high 
initial  momentum) 


Figure  66  Flow  Type  A3  -  Right  Coils  and  Final  Vortex 
(  Re  =  508,  V  =  150* ,  L'=  0,  DP-Lr-high 
initial  momentum  ) 
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Figure  67  Flow  Type  A2  -  Lower  Left  Coil  and  Final  Vortex 
(  Re  =  508,  V  =  120°  ,  L* =  0,  DP-Ll-high  initial 
momentum  ) 


Figure  68  Flow  Type  A3  -  Left  Coils  and  Final  Vortex 

(  Re  =  400,  V  =  1 55°,  L ' =  0,  DP-Ll-high 
initial  momentum  ) 


Figure  69  Flow  Type  AB  -  Lower  Coils 

(  Re  =  422,  y  =  90°  ,  L '  =  3. 15,  DP-Ki- 
upper  fig.  high  initial  momentum 
lower  fig.  low  initial  momentum  ) 


Figure  70 


Flow  Type  AB  -  Lower  Right  Coil 
(  Re  =  414,  V  =  90°  ,  L ’ =  3.15,  DP-Lr  ) 


107 


Figure  71  Flow  Type  AB  -  Lower  Left  Coil,  Two  Coil  Structu¬ 
re  in  the  Spacer  and  Final  Vortex 
(  Re  =  414,  t  =  90c  ,  L'  =  3-15,  DP-L1  and  AO) 
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Figure  72 


Flow  Type  B2  -  Large  V 
(  Re  =  422,  r  =  4 5°,  L'=  3*15,  DP-AQ) 
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Figure  73  Flow  Type  B2  -  Lower  Right  Coil  and  Final  Vortex 
(  Re  =  414,  4*  =  90°  ,  L’  =  3.I5,  DP-Lr  ) 
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Figure  74  Flow  Type  B2  -  Lower  Left  Coil 

(  Re  =  414,  V  =  90°  ,  L'  =  3.15,  DP-L1  and  AO) 
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Figure  75  Flows  Type  B2  -  Different  V 

(  left  fig.  Re  =  485,  V  =  45°,  L'=  3.15,DP-AO) 
(  right .  fig. Re  =  340,  ^=135°,  L'=  3.15,DP-H0C- 
high  initial  momentum  ) 
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Figure  ?6  Flow  Type  B2  -  Lower  Left  Coil  and  Final  Vortex 
(  Re  =  34o,  V  =  1350,  L’=  3.15,  DP-K1  ) 
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Figure  77  Flow  Type  B2  -  Lower  Right  Coil 

(  Re  =  3^0,  V  =  135°,  L’ =  3-15,  DP-Kr  ) 
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Figure  78  Final  Vortex 

(  Re  =  450,  DP-0  ) 


CHAPTER  VII 
CONCLUSIONS 

7  » 1  Summary 

In  conclusion,  for  the  investigated  range  of  variables, 
the  experimental  findings  can  be  summarized  as  follows: 

1.  lor  L  configuration,  there  exists  a  two  vortex  struc¬ 
ture,  symmetrical  with  respect  to  the  plane  of  curvature  of 
the  elbows. 

2.  i4or  ohe  S"  configuration,  there  exists  the  four  vortex 
structure  (for  Re  ^450),  symmetrical  with  respect  to  the 
plane  of  curvature  of  the  elbows. 

3*  4  or  the  angles  0  ^  V  (.  180  there  exist  three  basic  types 

Oj.  flow ,  identified  on  the  basis  of  the  two  coil  structure 
at  the  inlet  to  the  second  elbow. 

4.  For  the  three-dimensional  configuration,  there  is  a 
single  vortex  in  the  outlet  pipe  at  points  sufficiently 
downstream  from  the  second  elbow. 

5 •  There  exists  only  one  (second)  separation  region  when 
both  L  and  are  small,  and  there  exist  two  independent 
separation  regions,  one  in  each  elbow,  for  all  other  cases. 

6.  There  exists  backflow  inside  the  first  separation  zone 
and  for  the  higher  Reynolds  numbers  also  in  the  second  one. 

?.  The  backflow  in  the  first  elbow  is  amplified  by  the 
centrifugal  force  in  the  second  elbow  (especially  for 
*  =  180  )  . 

8.  The  backflow  in  the  second  elbow  changes  the  flow  pat¬ 
terns  downstream. 
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2- '.1 _ Suggestions  for  Further  Work 


As  has  been  previously  explained,  the  apparatus  has 
been  found  inadequate  to  give  valuable  quantitative  as  well 
as  qualitative  pitot  tube  confirmation  of  the  flow  patterns 
investigation  results  obtained  in  this  work.  Eecause  it 
is  doubtful  if  such  confirmation  can  be  obtained  for  the 
low  range  of  Reynolds  numbers,  the  more  advanced  measure¬ 
ment  techniques  are  required.  Eoth  the  laser  and  the  hot 
wire  anemometry  seem  to  be  suitable,  but  are  very  laborious 
and  require  the  specially  designed  apparatus  and/or  measur¬ 
ing  devices  (Nerem  et  al  1972,  Agrawal  et  al  1978)  . 

for  better  modelling  of  the  human  aortic  arch,  is  is 
suggested  to  change  the  shape  of  the  inlet  to  the  first  el¬ 
bow  in  order  to  avoid  the  separation  at  the  sharp  edge. 

Such  re-designed  apparatus  would  allow  to  investigate  the 
transition  to  turbulence  using,  for  example  a  hot  wire 
anemometer . 

The  experiments  using  the  dye  technique  should  be 
carried  out  with  a  larger  apparatus  in  order  to  make  the 
observations  easier  and  to  obtain  more  precise  quantitative 
results,  but  if  city  water  is  used,  also  an  efficient  me¬ 
thod  of  air  bubbles  removal  should  be  developed  and 
adequate  warm  water  source  should  be  provided. 

For  better  control  of  the  low-Reynolds  number  flow 
(  below  400  ) ,  it  is  suggested  to  use  a  precise  by-pass  valve 
at  the  inlet  to  the  inlet  tank. 
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TABLES  OF  RESULTS  FOR 
SECTION  5-2 
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Table  A.l  Dependance  of  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  0 


V  = 

0° 

V  = 

450 

<r= 

90° 

Re 

L^/mm/ 

Re 

L^/mm/ 

Re 

L^/mm/ 

j 

2637 

-4  .0 

1835 

-3.2 

from 

from 

140? 

0 

1206 

+  3-2 

2517 

not 

2637 

not 

1082 

+ 

• 

O 

to 

obs . 

to 

obs . 

879 

+  3*2 

450 

450 

645 

+4 . 0 

527 

]  not 

1 

i 

1  i 

1 

1 

1  1 

450 

J  obs . 

I 

1 

r  = 

135° 

y 

=  180° 

Re 

'  ] 

L1/mm/ 1 

1- 

Re 

.  ] 

41/mm/ 

2637 

-8.7 

2559 

-12.7 

1835 

-6.8 

1812 

-12.7 

1407 

-4.8 

1450 

-  9.5 

1206 

-2.3 

1242 

-9.5 

1082 

-1.6 

1012 

-7.9 

879 

-3.2 

870 

-7.1 

649 

-2.8 

544 

-9.5 

527 

-2.3 

444 

-9.5 

450 

1 

not  obs. 
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Table  A. 2  Dependance  of  L1  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L*  =  0.85 


Y  =  o*  V 


Re 

L1  (mm ) 

Re 

2675 

-9.5 

2482 

19^5 

-7.9 

1835 

142? 

-6.3 

1407 

1223 

-3.2 

1206 

1097 

-3.2 

1082 

856 

-2.4 

938 

594 

-5.6 

812 

428 

-5.6 

603 

496 

45* 

f  = 

90° 

L1  (mm) 

Re 

(mi 

-7.9 

2482 

-10.3 

-6.3 

1835 

-8.7 

-4.8 

1407 

-6.3 

-6.3 

1206 

-5.6 

-4.0 

1082 

-4.8 

-4.0 

981 

-4.8 

-4.8 

796 

-4.0 

-4.0 

603 

-3.2 

-8.7 

422 

-5.6 

V  = 

135° 

r  = 

180° 

Re 

L^(  mm) 

Re 

(mm) 

2482 

-10.3 

2562 

-11.9 

1835 

-10.3 

1208 

-9.7 

1407 

-4.8 

1367 

-8.7 

1206 

-4.0 

1139 

-8.7 

1082 

-4.0 

932 

-7.9 

938 

-8.1 

745 

-4.8 

812 

-6.3 

466 

-6.8 

603 

-2.4 

450 

-6.3 

469 

-1.6 

L 
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Table  A.  3  Dependance  of  L1  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  1.8 


*  = 

0° 

*  = 

45® 

V  = 

90° 

Re 

(mm ) 

Re 

(mm ) 

Re 

(mm ) 

2140 

-9.5 

2482 

-7.1 

2637 

-11.1 

1712 

-7.9 

1758 

-6.8 

1835 

-9.5 

1427 

-8.7 

1407 

-7.1 

1455 

-9.5 

1157 

-6.3 

1241 

-5.6 

1241 

-4.8 

951 

-4.8 

1055 

-3.2 

1111 

-2.1 

608 

-4.8 

861 

-1.6 

898 

-4.8 

535 

-4.8 

740 

0 

812 

-3.2 

441 

-4.0 

527 

-1 . 6 

603 

0 

457 

-1.6 

444 

-2.5 

V  = 

135° 

V 

=  180° 

Re 

(mm  ) 

Re 

(mm ) 

2637 

-9.5 

2718 

-11.1 

1835 

-8.7 

1812 

-7.9 

1407 

-7.1 

1450 

-7.9 

1241 

-4.8 

1176 

-7.1 

1111 

-6.3 

1087 

-7.9 

959 

-4.0 

763 

-11.1 

754 

-4.8 

512 

-9.8 

563 

-4.8 

473 

-6.3 

435 

-2.4 
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Table  A. 4  Dependence  of  ^  on  V  and  Re  for  a  Uniform  Velo- 
city  Profile  and  L’  =  3.15 


y 

=  0* 

f  = 

45. 

. .  -  _ 

v  = 

90° 

Re 

(mm  ) 

Re 

(mm  ) 

Re 

(mm ) 

2482 

j 

-11.1 

2482 

i 

-7.9 

1 

2675 

-10.3 

1758 

-9.5 

1835 

-7.9 

1861 

-9.5 

1407 

-7.9 

1407 

-7.1 

1585 

-8.7 

1206 

-6.3 

1241 

-6.3 

1337 

-7.5 

1110 

-6.3 

1141 

-5.0 

1157 

-5.6 

938 

-4.8 

938 

-3.2 

1070 

-4.8 

812 

-3.2 

639 

-2.4 

764 

-6.3 

496 

-3.2 

586 

-2.4 

542 

-5.6 

206 

-3.2 

449 

0 


460 

-4.8 

V  = 

135“ 

V  = 

180® 

Re 

(mm  ) 

Re 

Lj,  (mm) 

2517 

-11.1 

2559 

-11.1 

1783 

-10.3 

1812 

-7.9 

1427 

-7.9 

1450 

-7.9 

1223 

-6.3 

1242 

-7.9 

1070 

-4.8 

1061 

-6.3 

911 

-3.2 

791 

-4.8 

793 

-6.3 

473 

-7.9 

535 

-3.2 

476 

-3.2 
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Table  A. 5  Dependance  of  on  Y  and  Re  for  a  Uniform 

Velocity  Profile  and  L'  =  0 


Y  = 

0° 

1 

v  =  45*  I 

1 

1 

Y  = 

90° 

Re 

l 

L^/mm/ 

Re  Ly'mm/ 

1 

1  Re 

L^/mm/ 

2637 

+  9.5 

1835 

+9  •  5 

from 

from 

1407 

+8 . 2 

2517 

not 

2637  not 

1206 

+7  .8 

to 

obs . 

to  obs. 

1082 

+6.9 

450 

450 

879 

+6.3 

649 

+  5*6 

527 

]  not 

1 

450 

'  obs . 

1 _ 

I 

•  1 

1  ... _ 

.... _ 

_ t 

Y  = 

135* 

Re 

Ly'mm/ 

2637 

+  12.7 

1835 

+8.7 

1407 

+7.9 

1190 

+  7.9 

1082 

+6.9 

879 

+  6.3 

649 

+3*6 

527 

+5.6 

450 

pot  Obs.: 

Y  = 

►-* 

00 

0 

e 

i. 

Re 

L  q/ mm/ 

2559 

+  19.0 

1812 

+  22.2 

1450 

+19.9 

1242 

+  22.2 

1012 

+  20.8 

870 

+  12.7 

544 

+  12.7 

444 

_ 

not  obs 

127 


Table  A. 6  Dependance  of  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  0.85 


V  = 

0* 

T 

= 

45. 

- ■> 

V 

=  90* 

Re 

L^/mm/ 

Re 

L~/mm/ 

Re 

L^/mm/ 

-7.9 

2675 

-3.2 

2482 

-5. 

1 

1 

3 

2482 

1945 

-3.2 

1835 

-0  . 

5 

1835 

-7.1 

1427 

-7.1 

1407 

-4. 

8 

1407 

-7.1 

1223 

0 

1206 

+  9. 

5 

1206 

0 

1097 

0 

1082 

+4. 

3 

1082 

+  9.0 

856 

-1.7 

938 

+  1.2 

981 

+4.6 

594 

-4.2 

812 

-3. 

6 

796 

0 

428 

-6.4 

496 

-7. 

9 

603 

-3.9 

i - 

1 

V 

.  _  - 1 

422 

i 

1 .. 

-7.9 

V 

=  135 

0 

r 

! 

| 

f  = 

180# 

Re 

L^/mm/ 

! 

Re 

L^/mm/ 

2482 

+  9. 

5 

2562 

+6 . 3 

1835 

+9. 

5 

1708 

+6 . 3 

1407 

+9. 

5 

1367 

+6 . 3 

1206 

+9. 

5 

1139 

+6.3 

1082 

+  12 

.1 

932 

+7.9 

938 

+  14 

•  3 

745 

+7.9 

812 

+  11 

.1 

466 

+7.9 

603 

+  5  • 

3 

450 

+6 . 3 

469 

-7. 

.  i 

9 

-  - 

128 


Table  A. 

7  Dependanc 

e  of  Lj 

on  Y  and  Re  for 

a  Uniform 

Velocity 

Profile 

and  L'  =1.8 

V 

=  0° 

i 

f 

0 

vn 

-3- 

11 

V 

=  90" 

Re 

( 

L^/mm/ 

Re 

Lj/mm/ 

Re 

Lj/mm/ 

2140 

-28.6 

2482 

-28.6 

2637 

-28.6 

1712 

-26.9 

1758 

-19.0 

1835 

-25.4 

1427 

-23.4 

1407 

-19.0 

1455 

-22.2 

115? 

-21.7 

1241 

-19.0 

1241 

-19.5 

951 

-19.1 

1055 

-19.0 

1111 

-9.5 

808 

-23.8 

861 

-24.2 

898 

-12.5 

535 

-28.6 

740 

-28.6 

812 

-16.1 

441 

-30.2 

527 

-36.5 

603 

-19.5 

1 

457 

-36.5 

444 

-21.2 

Y  = 

135° 

0 

0 

CO 

1 — 1 

II 

Re 

Ly'mm/ 

Re 

L^/mm/ 

2637 

-28.6 

2712 

-28.6 

1835 

-24.1 

1812 

-25.3 

1407 

-19.5 

1450 

-22.7 

1241 

-9.5 

1176 

-19.1 

1111 

-9.5 

1087 

-9.5 

959 

-19.5 

763 

+  6.3 

754 

-19.5 

512 

+6.3 

563 

-22.2 

473 

+6 . 3 

435 

-22.2 

l 

1 _ 

i  _ _ _ _ _ 

129 


Table  A. 8  Dependence  of  on  ¥  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  3.15 


¥ 

=  0° 

¥  = 

•£- 

V_n. 

0 

¥  = 

90® 

Re 

L^/mm/ 

Re 

L^/rnm/ 

Re 

L^/mm/ 

2482 

-73.2 

2482 

-7  4.6 

2675 

-58.3 

1758 

-73.2 

1835 

-70  .6 

1861 

-59.4 

1407 

-730 

140? 

-67.5 

1585 

-60.3 

1206 

-73.2 

1241 

-65 . 1 

1337 

-57-6 

1110 

-73.2 

1141 

-62.8 

1157 

-60.3 

938 

-68  .7 

938 

-60.3 

1070 

-60.3 

812 

-60  .4 

639 

-62.8 

764 

-60.3 

496 

-60  .4 

586 

-65.I 

542 

-60.3 

206 

-47.7 

449 

not  obs. 

460 

-58.3 

!' 


¥  = 

135° 

¥  = 

180° 

Re 

Ly4nm/ 

Re 

L~/ mm/ 

2517 

-66 .7 

2559 

-74.6 

1783 

-63.4 

1812 

-68.4 

1427 

-60.3 

1450 

-60.3 

1223 

-57.0 

1242 

-60.3 

1070 

-53.5 

1061 

-60.3 

911 

-57.0 

791 

-47.7 

793 

-58.3 

473 

-47.7 

535 

-60.3 

4  76 

-57.7 
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Table  A. 9  Dependance  of  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  0 


V  = 

0° 

Y  = 

j _ _ _ 

450 

f  = 

90* 

Re 

L  £./ mm/ 

Re 

L^/mm/ 

Re 

LVmm/ 

2517 

-5.6 

2637 

-6.3 

2637 

-6.3 

1783 

-3.2 

1835 

-4.8 

I835 

-5.2 

1476 

0 

1407 

-3.2 

1407 

-3.2 

1259 

0 

1206 

0 

1206 

-1.6 

1157 

0 

959 

+6.3 

1082 

0 

839 

+  11.1 

767  ■ 

not 

879 

0 

522  ' 

not 

544 

* 

649 

+6 . 3 

450 

obs . 

450 

obs . 

527  1 

not 

450  J 

obs . 

=  1 35“ 

Re 

L^/mm/ 

2637 

-10.3 

1835 

-  9.5 

1407 

-5.6 

1206 

-4.8 

1082 

-3.2  ! 

879 

0 

649 

+  3.2 

527 

+9.5 

450 

not  obs. 

j 

[ . r: 

180° 

Re 

L^/ mm/ 

2559 

-15.9 

1812 

-15.3 

1450 

-15.9 

1242 

-15.9 

1012 

-14.3 

870 

-15.9 

544 

-15.9 

444 

; 

1 

| 

-11.1 
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Table  A. 10  Dependence  of  on  V  and  Re  for  a  Uniform 

Velocity  Profile  and  L'  =  0.8 5 


V  = 

0° 

V 

Re 

L^/mm/ 

Re 

26?5 

-4.8 

2482 

1945 

-4.8 

I835 

142? 

0 

140? 

1223 

+  2.4 

1206 

1097 

+4.8 

1082 

856  ' 

not 

938 

594 

» 

obs . 

812 

428  J 

603 

I 

i  496 

i 

|  ! 

i.  t 


135° 

Re 

L^/mm/ 

2482 

-13.5 

1835 

-9.5 

1407 

-7.9 

1206 

-10.3 

1082 

-7.1 

00 

ON 

-6.3 

812 

”5  •  6 

603 

+3.2 

469 

lot  obs. 

i 

i 


45° 

L./mm/ 

V  = 

Re 

90* 

L^/mm/ 

-7.9 

2482 

-11.1 

-6.3 

1835 

-7.9 

-1.6 

1407 

-4.8 

0 

1206 

-2.4 

0 

1082 

-1.6 

+3.2 

981 

-2.4 

+9.5 

796 

+4.8 

not 

603 

not 

obs . 

422 

obs . 

^  = 

=  180® 

Re 

L  mm/ 

2562 

-15.9 

1708 

-15.9 

1367 

-15.9 

1139 

-15.9 

932 

-13.5 

745 

-15.9 

466 

-6.3 

450 

-4.8 

i 
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Table  A. 11 


Dependance  of  on  V  and  Re  for  a  Uniform 
Velocity  Profile  and  L'  =  1.8 


V  = 

0* 

V  = 

45° 

V 

=  90° 

Re 

L^/mm/ 

Re 

L^/mm/ 

Re 

L  r/ mm/ 

2140 

-4.8 

2482 

-8.3 

2637 

-9.5 

1712 

0 

1758 

-6.3 

1835 

-8.3 

1427 

0 

1407 

-2.4 

1455 

-7.1 

115? 

+  1 .6 

1241 

-2.4 

1241 

-4.0 

951 

+4 . 0 

1055 

0 

1111 

0 

808 

not 

861 

+6 . 3 

898 

0 

535 

* 

obs . 

740 

not 

812 

+6.3 

441 

527 

obs . 

603 

+9.5 

j _ 

457  J 

1 . .. 

444 

not  obs 

V 

=  135# 

V 

=  180° 

Re 

L^/mm/ 

Re 

L^/mm/ 

2637 

-9.7 

2718 

-12.7 

1835 

-9.5 

1812 

-11.1 

1407 

-7.9 

1450 

-9.5 

1241 

-7.1 

1176 

-14.3 

1111 

-4.0 

1087 

-13.3 

959 

-6.3 

763 

-11.1 

764 

-1.6 

512 

-6.3 

603 

+9.5 

473 

0 

435 

not  obs. 

i 
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Table  A. 12  Dependance  of  L5  on  V  and  Re  for  a  Uniform 

Velocity  Profile  and  L*  =  3. 15 


V 

=  0° 

v  = 

450 

« 

0 

o\ 

11 

Re 

Ly'rnm/ 

Re 

L  £./ mm/ 

Re 

L^/ mm/ 

2482 

-7.9 

2482 

-90 

2675 

-9.5 

1758 

-6.3 

1835 

-7.8 

1861 

-7.9 

1407 

-4.0 

1407 

-5.6 

1585 

-6.3 

1206 

0 

1241 

-3.2 

1337 

-5.6 

1110 

0 

1141 

-1.6 

1157 

-4.0 

938 

0 

938 

-1.6 

1070 

-2.4 

812 

+5.0 

639 

+7.9 

764 

-0  .8 

496 

+11.1 

584 

not 

542 

+6.3 

206 

not  obs . 

449 

r 

obs . 

46o 

not  obs . 

V  = 

135“ 

V 

=  180° 

Re 

L^/mm/ 

Re 

L^/mm/ 

2517 

-10.3 

2559 

-14.3 

1783 

-5.6 

1812 

-8.9 

1427 

-9.5 

1450 

-10.3 

1223 

-4.8 

1242 

-8 .5 

1070 

-4.0 

1061 

-6.3 

911 

-4.8 

791 

-9.5 

793 

-1.6 

473 

0 

535 

+4.0 

476 

not  obs. 
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Table  A.13  Dependance  of  L1  on  L',  V  and  Re  for  a  Fully 

Developed  Velocity  Profile 


- 

11  1 

.° 

V  =  15* 

L'=  0,  V  =  90° 

- - 

L'  =  0, 

Y  =  180* 

Re 

L^/mm/ 

Re 

L^/mm/ 

Re 

L^/ mm/ 

2500 

not 

- - 

2500 

-19.0 

1835 

obs . 

1835 

-19.0 

1450 

1450 

-15.9 

1240 

+  10.3 

1240 

-14.3 

from 

1110 

+7.9 

1110 

-11.9 

2500 

not 

900 

+  6.3 

900 

0 

to 

obs . 

812 

+9  •  5 

812 

0 

1*53 

650 

+12.7 

680 

-1 .6 

*- 

_ 

' 

453  not  obs. 

48  5 

0 

L'=3.15 

,  V  =180® 

L'=3-l 5 

V  =90° 

1 L’=3.15 

,  =15c 

Re 

L^/mm/ 

Re 

L^/mm/ 

Re 

L^/mm/ 

2500 

-19.0 

2500 

-15.9 

2500 

-4.0 

1835 

-15.9 

1835 

-9.5 

1835 

-3.2 

1450 

-11.1 

1450 

-4.8 

1450 

-3.2 

1240 

-9.5 

1240 

-5.6 

1240 

-3.2 

1110 

-7.9 

1110 

-4.8 

1110 

-5.6 

900 

-7.1 

900 

-4.8 

900 

-4.0 

812 

-4.0 

812 

-1.6 

812 

-4.8 

650 

-2.4 

650 

-0.8 

550 

-1.6 

612 

0 

453 

-4.8 

444 

.  i 

-3.2 

135 


Table  A. 14  Dependance  of  on  L',  V  and  Re  for  a  Fully 

Developed  Velocity  Profile 


L ' =  0,  V  =  150 


r  “  "» 

Re 

L^/iran/ 

from 

2500 

not 

to 

obs . 

CO 

-V 

L 

;L’=3.i5, 

V  =15° 

Re 

L^/mm/ 

2500 

-68.2 

1835 

-64.3 

1450 

-64.3 

1240 

-74.6 

i  1110 

-68.2 

900 

-74.6 

812 

-74.6 

650 

-68.2 

453 

-68.3 

L'=  0,  * 

=  90° 

Re 

L^/mm/ 

2500 

not 

1835 

► 

obs . 

1450  . 

1240 

+  12.7 

1110 

+  12.7 

900 

+  12.7 

812 

+  12.7 

650 

+  12.7 

453 

not  obs. 

L’=3.15 

,  V  =90c 

Re 

L^/mm/ 

j 

2500 

-73.2 

1835 

-7  4.6 

1450 

-74.6 

1240 

-74.6 

1110 

-74.6 

900 

-74.6 

812 

-74.6 

650 

-68.2 

453 

-68.2 

L  '  =  0, 

V  =  180° 

Re 

Ly'mm/ 

2500 

+  25.4 

1835 

+  19.0 

1450 

+  14 .9 

1240 

+  19.0 

1110 

+  20.6 

900 

+  12.7 

812 

+  15.9 

680 

+  12.7 

48  5 

+  12.7 

j  L •  =  3 . 15 ,  V  =180° 


Re 

Ly6nm/ 

2500 

-68.3 

1835 

-74.6 

1450 

-74.6 

1240 

-7  4.6 

1110 

-74.6 

900 

-74.6 

812 

-74.6 

550 

-68.2 

444 

-60.3 

1 

_ J. 
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Table  A. 15  Dependance  of  on  L',  Y  and  Re  for  a  Fully 

Developed  Velocity  Profile 


L' =  0,  V  =  15° 


Re 

L^/mm/ 

2500 

-7.9 

1835 

-7-1  | 

1450 

-3.2 

1240 

0 

1110 

+19.0  ! 

900 

812 

not 

650 

► 

obs . 

453 

L'  =  0, 

=  90° 

Re 

1 

L^/mm/ 

j 

2500 

• 

!  -4 

1835 

-3.2 

1450 

0 

1240 

+2.4 

1110 

+2.4 

900 

+6 . 3 

812 

+9.5 

650 

1  not 

453 

j  obs . 

i 


'  =  0, 

V  =  180* 

Re 

L^/mm/ 

2500 

-14.3 

1835 

-14.3 

1450 

-15.9 

1240 

-15.9 

1110 

+  14.3 

900 

+  13.5 

812 

-11.1 

680 

-9.5 

485 

-4.8 

L'  =  3-15 

V  =150 

Re 

] 

L^/mra/ 

j 

2500 

j  1835 

-7.1 

1450 

-4.0 

j  1240 

-1.6 

1110 

0 

900 

0 

812 

+4 . 8 

650 

I  not 

453 

1  obs . 

L’  =  3.15 

,  Y  =90° 

Re 

L^/mm/ 

2500 

-12.7 

1835 

-8.0 

1450 

-4.0 

1240 

-2.4 

1110 

-0.8 

900 

-1.6 

812 

-1.6 

650 

+7.9 

453 

- 

not  obs. 

L,=3-15 

1 

,  ^  180° 

Re 

L^/mm/ 

2500 

-15.9 

1835 

-8.5 

1450 

-4.0 

1240 

0 

1110 

0 

900 

-1.6 

812 

0 

550 

0 

444 

+  2.4 

